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1.  SUMMARY 


This  technical  report  summarizes  the  theory  based  mathematical  models,  describing  various 
processes  that  occur  in  an  operating  lithium-based  cell.  Such  models  provide  the  guidelines  for 
cell/battery  design  improvement;  leading  to  the  experimental  work  reduction  to  attain  a  desired 
cell  design  for  optimum  cell  performance.  The  developed  formulation  covers  the  following  sub 
topics: 

1.1  Cathode  Formulation 

The  developed  mathematical  expressions  describe  the  transient  intercalated  lithium  concentration 
profiles  in  the  cathode  particles  of  an  active  material,  for  example,  FePc  (iron  phthalocyanine), 
transient  cathode  current  and  charge  storage  during  the  cell  discharge  period.  Voltage  loss 
associated  with  the  lithium  ion  and  electron  transport  in  the  cathode  composite  of  electrolyte, 
carbon  black  and  active  material  can  be  estimated.  Based  on  the  thermodynamic  equilibrium 
condition  for  the  electrochemical  reaction  at  the  electrolyte-active  material  interface,  a 
mathematical  equation  relates  the  cathode  reaction  reversible  electric  potential  to  the  lithium 
content  of  the  cathode  electrode.  Based  on  the  fundamentals  of  electrode  kinetics,  the  cathode 
active  material-electrolyte  interfacial  current  density  has  been  quantitatively  related  to  the 
cathode  active  material  lithium  content,  lithium  ion  concentration  in  the  electrolyte,  the  surface 
overpotential  and  the  cathodic  reaction  exchange  current  density. 

1.2  Electrolyte  Separator  Formulation 

To  compute  the  voltage  drop  in  the  separator,  e.g.,  LiBF4  in  microporous  poly(ethylene  oxide), 
between  the  cell  electrodes,  a  relatively  simple  mathematical  expression  is  provided. 

1.3  Anode  Formulation 

The  electrokinetics  of  the  lithium  metal  sheet  in  contact  with  the  electrolyte  has  been  modeled  to 
predict  the  interfacial  current  density.  The  interfacial  current  density  is  expressed  in  terms  of  the 
anode  exchange  current  density  and  the  surface  overpotential  at  a  given  temperature  at  the 
anode-electrolyte  interface.  A  developed  mathematical  expression  relates  the  reversible  electric 
potential  of  the  anode  reaction  to  its  standard-state  electric  potential  and  the  lithium  ion  activity 
in  the  electrolyte  at  the  lithium  anode-electrolyte  interface. 

1.4  The  Lithium-Based  Cell  Reversible  Electric  Potential  Formulation 

The  developed  formulation  relates  the  reversible  cell  potential  to  the  lithium  activity  in  the 
cathode  active  material,  e.g.,  lithium  in  MyPc  (e.g.,  FePc).  Using  the  developed  mathematical 
expression,  one  can  develop  semiempirical  relation  between  the  lithium  activity  in  the  cathode 
active  material  and  lithium  content  of  the  cell  cathode  from  the  experimental  data  on  the  cell 
reversible  voltage  versus  the  cell  lithium  content.  Such  an  empirical  relation,  in  turn,  may  be 
used  to  predict  the  actual  cell  voltage  during  the  cell  discharge  period. 
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1.5  Cell  Performance  Formulation 


The  actual  cell  voltage,  electric  power  and  thermodynamic  efficiency  are  predicted  during  a 
given  cell  discharge  period. 

1.6  Cathode  Active  Material  Spherical  Particle  Utilization  to  Store  Lithium  during  the 
Cell  Discharge  Period 

The  developed  formulation  predicts  the  transient  lithium  concentration  profile  in  a  spherical 
particle  of  the  cathode  active  material.  The  formulation  was  developed  for  a  constant  active 
material-electrolyte  interfacial  current  density  with  the  assumption  of  lithium  diffusion  inside  a 
spherical  particle  of  the  cathode  active  material  controlling  the  overall  cathode  and,  hence,  the 
overall  cell  performance  during  the  cell  discharge  period.  Furthermore,  a  developed 
mathematical  expression  predicts  the  cathode  active  material  utilization  or  its  coulomb  efficiency 
until  the  lithium  concentration  in  the  cathode  active  material  phase  at  its  interface  with 
electrolyte  phase  reaches  the  saturation  level  when  the  actual  cell  voltage  rapidly  decreases. 

1.7  Voltage  Loss  associated  with  the  Diffusion  Transport  of  Lithium  in  the  Cell  Cathode 
Active  Material 

The  Nemst-Einstein  relation  linking  the  lithium-ion  mass  diffusivity  and  its  ionic  conductivity 
was  employed  to  develop  a  mathematical  expression  for  the  prediction  of  transient  voltage  loss 
associated  with  the  lithium  diffusion  in  the  cell  cathode  active  material. 

1.8  Effect  of  Lithium  Storage  on  the  Cell  Reversible  Voltage 

This  sub  section  briefly  describes  the  relation  showing  the  transient  effect  of  lithium  storage  in 
the  cathode  active  material  on  the  cell  reversible  voltage  during  the  cell  discharge  period. 

1.9  Further  Treatment  of  the  Reversible  Potential  of  the  Lithium-Based  Cell 

The  liquid-phase  electrolytic  solution  theory  has  been  extended  to  account  for  the  nonideal 
behavior  of  the  solid-state  cathode  active  material  containing  the  intercalated  lithium  in  it 
through  the  relation  between  the  lithium  ion  content  of  the  cathode  active  material  and  lithium 
activity  coefficient.  This  information,  in  turn,  has  been  related  to  the  cell  reversible  voltage. 

1.10  Fractional  Utilization  of  a  Cathode  Active-Material  Particle  as  a  Function  of  Time 

In  this  brief  sub  section,  the  formulation  has  been  provided  to  compute  the  fractional  utilization 
of  the  cathode  active  material  in  a  spherical  particle  as  a  function  of  time  at  a  constant  active 
material-electrolyte  interfacial  current  density  during  the  cell  discharge  period. 
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2.  INTRODUCTION 


The  lithium-based  battery  development  is  being  pursued  for  aerospace,  defense,  medical,  and 
hybrid  electric  and  battery  electric  vehicle  applications.  Lithium-based  storage  devices  also 
provide  power  for  laptop  computers,  cellular  phones,  camcorders,  etc.  In  the  excellent  review  on 
lithium  batteries  and  cathode  materials  [1],  it  has  been  stated  that  the  layered  oxides,  e.g., 
LixNio.4Mno.4Coo.2O2,  continued  their  current  improvement  with  mixed  transition  metals  slowly 
replacing  the  pure  cobalt  system  (i.e.  LiCo02)  as  the  cathode  materials  for  the  lithium-based 
batteries.  They  are  expected  to  offer  enhanced  safety  with  enhanced  capacity  at  a  lower  cost.  The 
first  potentially  low-cost  cathode  material  appears  to  be  the  lithium  iron  phosphate  (LiFeP04). 
This  will  be  achieved  only  if  the  costs  of  the  other  cell  components  are  also  reduced.  These 
include  the  electrolyte,  a  lower  cost  anode  with  enhanced  safety  features  that  include  the  use  of 
graphite  and  the  use  of  ‘thicker  cathodes’  to  reduce  the  cost  of  electrode  support  materials, 
separators,  etc.  Metal  phthalocyanines  (e.g.,  Li2Pc[LLC32Hi6N8],  [2])  are  also  being  investigated 
as  low  cost  cathode  materials. 

Cost,  specific  energy,  power,  cycle  life,  calendar  life,  and  safety  are  important  factors  that  must 
be  considered  while  deciding  a  lithium-based  system  for  a  particular  application.  Selection  of 
component  materials  would  decide  cost  and  life  of  a  lithium-based  system.  The  required  specific 
energy  and  power  of  a  battery  depend  on  the  cell  design.  Experimental  determination  of  an 
optimum  cell  design  to  provide  a  desired  cell  performance  requires  an  intensive  work  that 
involves  fabricating  lithium-based  cells  of  various  designs,  e.g.,  various  insertion  electrode 
compositions  (such  as  the  intercalated  lithium  metal  oxide,  lithium  metal,  intercalated  lithium-in 
carbon,  etc.),  their  thicknesses,  void  fractions,  electrode  active  material  particle  sizes,  an 
electrolyte  composition,  and  thickness  of  an  electrolyte-containing  separator  between  the  cell 
electrodes.  Theory  based  mathematical  models,  describing  various  processes  which  occur  in  an 
operating  lithium-based  cell,  are  very  helpful  in  providing  the  guidelines  for  cell/battery  design 
improvement;  leading  to  the  experimental  work  reduction  to  attain  a  desired  cell  design. 

A  typical  lithium-based  cell,  lithium  I  electrolyte  LiBF4  in  poly(ethylene  oxide)  I  intercalated 
lithium  in  MyPc,  is  shown  in  Figure  1. 


Electrolyte 

Cathode 

Anode 

Separator 

Current 

(LI) 

(UBF4  in 

_ -  . 

Collector 

PEO) 

(Al) 

Figure  1.  Sketch  of  a  lithium-based  cell  (not  to  a  scale). 


The  various  cell  components,  such  as  the  anode,  electrolyte  separator,  cathode  and  cathode 
current  collector,  shown  in  Fig.  1,  are  explained  in  the  Fist  of  Acronyms,  Abbreviations,  and 
Symbols. 
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3.  METHODS,  ASSUMPTIONS,  AND  PROCEDURES 


The  mathematical  model  formulation  of  a  lithium-based  cell  developed  using  the  principles  of 
thermodynamics,  transport  phenomena,  and  electrochemical  kinetics  is  presented  in  the  next 
section. 
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4.  RESULTS  AND  DISCUSSION 


4.1  Mathematical  Model  Formulation 

4.1.1  Cathode  Formulation 

It  is  assumed  that  the  lithium  ion  transport  through  the  solid  electrolyte  to  the  active  material 
particles  of  metal  phthalocyanine  such  as  iron  phthalocyanine  (FePc)  or  dilithium  phthalocyanine 
is  much  faster  than  the  lithium  ion  diffusion  in  the  active  material  particles  of  the  cell  composite 
cathode.  A  lithium  ion  on  its  arrival  at  the  three  phase  interface  between  the  electrolyte,  metal 
phthalocyanine  (MyPc),  and  carbon  black  of  electronic  conductivity  much  higher  than  that  of 
MyPc  undergoes  a  fast  electrochemical  reaction:  Li++e~n  (Li=Li++e‘  in  the  cathode  active 
material).  Then,  lithium  (as  (Li++e  ))  is  assumed  to  enter  a  channel  between  two  layers  of  MyPc 
and  diffuses  in  a  solid-state  cathode  particle  of  MyPc  during  the  period  of  lithium  storage  in  the 
cathode  material,  i.e.  during  the  cell  discharge  period.  The  overall  process  of  lithium  storage  in 
the  cathode  active  material  is  assumed  to  be  controlled  by  the  lithium  diffusion  in  the  cathode 
MyPc  particles.  The  lithium  accumulation  is  assumed  to  take  place  under  the  transient, 
isothermal  and  isobaric  conditions.  The  differential  model  equation  describing  the  lithium 
diffusion  and  accumulation  in  a  spherical,  active  material  particle  is  represented  by  the  following 
derived  equation: 


dc  1  d  ( 
dt  r 2  dr 2  dr  , 


(1) 


where  c  =  molar  concentration  of  accumulating  or  intercalating  lithium,  D  =  molecular  mass 
diffusivity  of  lithium  in  the  active  material  particles.  The  initial  and  boundary  conditions  to  solve 
the  partial  differential  equation  (1)  are: 

The  initial  condition,  t  =  0;  c  =  a,  initial  uniform  concentration  of  lithium  in  the  cathode  active 
material  channels  in  a  spherical  particle  of  radius,  a. 

The  boundary  conditions  for  t  >  0: 

At  the  spherical  particle  center,  i.e.,  at  r  =  0,  c  has  a  finite  value;  and  at  the  particle  surface,  i.e. 
at  r  =  a,  c  =  c0. 

For  the  assumption  of  D  =  D  (composition  averaged  diffusivity)  =  constant,  the  partial 
differential  equation  (1)  becomes 


dc  -( d~c  2  dc 

—  =  D  — ^  + - 

dt  ydr~  r  dr  y 


(2) 


The  analytical  solution  of  the  differential  equation  (2)  using  the  above  mentioned  boundary 
conditions  [3]  leads  to: 
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The  transient  concentration  at  the  particle  center  can  be  obtained  from  Eq.  (3)  using  L’HopitaTs 
rule: 
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The  dimensionless  radial  distance,  E, ,  and  the  dimensionless  time,  r,  are  now  defined  as 
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The  equations  (3)  and  (4)  are  now  expressed  as 
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To  predict  the  stored  or  intercalated  lithium  concentration  in  a  spherical  particle  of  a  cathode 
active  material  for  short  times,  the  following  equation  may  be  used. 


(  c~ci  1 


\Co  C\  J 


erfc^tl-erfc^2n  +  1)+^ 


2  yfr 


(9) 


where 


erfc 


(2/7  +  1)-^ 


l- erf 


(2/7  +  1)-^ 
2  yfr 


(10) 
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erfc 


(2n  +  l)  +  <f; 
2  sfr 


=  \-erf 


(2n  +  l)  +  ^ 

i4t 


(11) 


where 


du 


eff{x)  =  ^=  J  e 

dx  u=o 
erfc  (0)  =  1,  erfc  (oo)  =  0 


(12) 

(13) 


Equations  (7)  and  (9)  can  be  used  to  predict 


f  \ 

c  —  c 

- —  as  a  function  of  ( r  and  d; ).  Equation  (8) 

vc*-cJ 


can  be  used  to  predict 


c-c, 


\Co~Cl  J 


as  a  function  of  r . 


f=o 


The  total  flux  of  lithium  into  a  spherical  particle  of  radius,  cl  at  any  time  is  given  by 


N1  =4  kcC 

1  v  panicle- flux  u 


D 


dc 

dr 


r=a  J 


(14) 


The  mathematical  expression  for 
(14)  to  obtain 


r dc 

\ 

Kdr 

r=a  J 

was  obtained  from  Eq.  (7)  and  substituted  into  Eq. 


N  particle— flux  =  [(8^aD)(c0-Cj)]5]exp(-nVr) 

n= 1 


(15) 


The  lithium  g-atom  balance  was  applied  over  a  spherical  particle  to  obtain  the  following 
expression  for  the  determination  of  stored  or  intercalated  lithium  at  any  time  during  the  cell 
discharge  period. 


N'=N' 


t=t 


|  dN'=(N'-Nl)=  \N'r.M^dt  (16) 

N'=N'0  t= 0 

The  use  of  Eq.  (15)  in  Eq.  (16)  leads  to: 

(/y'  -  iV')  =  (amount  of  lithium  stored  in  a  particle  of  the  cell  cathode  active  material  in  time, 
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(17) 


(jV'-AC): 


—  net'  (co  —  c, ) 


i-4z< 

n= 1 


exp(-n2;r2r)  j 
n2 


In  Eq.  (17),  N‘o  =  g-atom  (or  g-mole)  of  lithium  in  a  spherical  particle  of  the  cathode  active 
material  just  before  the  cell  discharge  to  store  or  intercalate  lithium  into  the  cell  cathode  and 
provide  electric  power  to  an  external  circuit.  The  total  amount  of  lithium,  N‘ ,  present  in  the 
active  material  particle  until  time,  r,  is  given  by 


N1  =N'  + 


-7ra\c0-Cl) 


i-4z< 

n= 1 


exp  |-nV2rj 


n 


(18) 


As  t —>  oo ,  therefore,  r  — >  oo ;  Eq.  (18)  reduces  to 


Nl  ~  N'0  +  —  kci3  (c0  -  Cj ) 


(19) 


For  the  situation  of  ‘no-storage’  of  lithium  in  the  cathode  active  material  channels  before  the  cell 
discharge,  the  initial  concentration  of  lithium,  c,  =  0  and  N’o  =  0  ;  then,  Eqs.  (18)  and  (19) 
reduce,  respectively,  to 


N'  = 


f4  3  ^ 

—  na  c 

\3  °J 


n= 1 


exp^-nW) 


n 


(20) 


Nl=—  nac„ 


(21) 


The  active  material  channel-stored  lithium  concentration,  c,  is  related  to 

g  -mole  of  the  channel-stored  or  intercalated  lithium 

x  = -  by  the  following  equation: 

g-mole  of  the  cathode  active  material  at  its  initial  state 


c  = 


f  Pm^Pc, 


.) 


V  Mm>Pc  j 


X 


(22) 


where  pu  Pc  0  =  density  of  a  cathode  active  material,  e.g.,  FePc  or  Li2Pc;  and  Mm  Pc  = 

molecular  weight  of  MyPc.  For  the  situation  of  cathode  active  material  expansion  during  the  cell 
discharge  period,  the  channel-stored  or  intercalated  lithium  concentration  c  and  x  are  related  by 
the  following  equation 


X 


(23) 


c  = 


Art  Pc,  c 


M 


MyPc  V 


X  3 
l  +  axy 


where  pM  Pc  0  =  density  of  the  cathode  active  material  when  the  channel-stored  lithium 

concentration  is  zero,  i.e.  x  =  0;  and  a  =  coefficient  of  the  active  material  particle  volume 
expansion. 

The  factor  (l  +  ax)  in  Eq.  (23)  is  related  to  the  active  material  density  by  the  following 
equation 

(l  +  ax)  AMyPc,  o 
l  AV[yPC 

where  pMPc  —  density  of  the  ‘initial’  active  material,  MyPc,  when  the  active  material 
channel-stored  lithium  content  is  x. 

If  the  initial  stored  lithium,  xini ,  were  known;  then,  the  initial  concentration  ,  a,  of  the 
intercalated  lithium  may  be  computed  from 


AiVlyPC,  O 


Mx 


v !  +  axini  j 


(25) 


If  it  is  assumed  that  the  stored  or  intercalated  lithium  in  a  cathode  active  material  particle  at  the 
particle-electrolyte  interface,  co  is  equal  to  cMd  ,  the  stored  lithium  concentration  when  the 

active  material  channels  are  full  of  lithium;  then  co  can  be  determined  from 


sat'd  Am  Pc,  o 

Ca  =  c  =  — - — 

M 

X/T  D,. 


sat'd  \ 


V1  +  ocxsa,'d  J 


(26) 


The  active  material,  lithium-diffusion  controlled  current  per  active  material  particle  of  radius,  a , 
is  given  by 


Ic  =  FN‘  .  (27) 

p  particle— flux  v  7 

Using  the  information  in  Eq.  (15)  for  N'particle_flux ,  Eq.  (27)  leads  to: 

n= co 

Ip  =F{SxaD(c0  -Cj)}  ^exp(-n2^2r),  [A]  (28) 

n= 1 

where  F  —  Faraday’s  constant  =  96ASlC/g-equiv. 
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The  cathodic  current  per  unit  external  surface  area  of  an  active  material  particle  of  initial  radius, 
V»’is  given  by 


A  2  F  A 


V  api’°  J 


/  2  2  \ 

r  a  i 

(-n  n  Tj , 

i 

5 

Vi  to 

1 _ 

(29) 


For  the  situation  of  constant  active  material  particle  volume  (i.e.  a  —  0),  Eq.  (29)  becomes 


•c 


f  2  F  A 


p,a= 0 


\api,oJ 


n=\ 


(  2  2  \ 

r  a  i 

(-n  71  T  J, 

i 

<N  03 

b 

_ 1 

(30) 


For  short  times,  the  cathodic  current  per  unit  external  surface  area  of  an  active  material  particle, 
iCp, shorn  is  given  by 


p, short 


FZ^^(l+q,pg 


pi,o 


n= 0 


erfc 


^  n  + 1^ 


yfr 


-erfc 


^  n  ^ 

\/r 


2A  ,  „-(»+1)' 


,2/,  \ 


sfnT 


cm. 


(31) 


The  diffusion-controlled  cathodic  current  per  unit  initial  mass  of  the  active  material  particles  of 
radius,  api  o ,  is  given  by 


or 


L  p,  per  unit  mass  ( ^  \ 


3^api,o 


Pm„Pc,o 


(32) 


Tc 

p,  per  unit  mass 


6  FD 

P MyPC,0  J 


(c0-cl), 


£exp(-nVr J, 
n= 1 


gm  of  active  material  particles  of  initial  size  a  . 


pi,o 


(33) 


for  the  assumption  of  a  —  0,  where  r  •  0  =  — —  . 

a 

pi,o 

For  the  situation  of  the  cathode  active  material  particle  volume  change,  Ic  perunitmass  per  unit 
initial  cathode  active  material  mass  may  be  estimated  from  the  following  equation: 


to 


(34) 


6FD 

1 

(co  -  c\ ),  (l  +  ax)3 

(  \ 
2  2 

~n  *  Tpi,o 

yPMyPc,o  j 

2  Z^ex  P 

api,o  n= 1 

2 

U l+ax)3  J 

The  total  transient  cathodic  current  associated  with  the  active  material  particles  of  size,  api  0 
present  in  the  entire  cathode  is  given  by 


l[c  = 


(  -  ,r\ 

GFDm'f 

( wPiAco-c \)A 

PmvPc,o 

K  api,o  J 

n= co 

Zexp  {-n27i2Tpio),[A\ 

n= 1 


(35) 


where  m'f  =  initial  amount  of  active  material  particles  of  all  sizes  in  the  cell  cathode  when  the 
channel  stored  lithium  amount  is  zero,  i.e.,  x  —  0\  and  wpj  0  =mass  fraction  of  active  material 
particles  of  radius,  api  0 ;  if  there  were  ‘no  intercalated  lithium’  in  the  cathode  active  material 
particles. 

If  the  presence  of  different  active  material  particle  size  ranges  is  accounted  for,  the  total 
lithium-diffusion  controlled  current  in  the  cell  cathode  is  given  by 


jt,C  _ 


(  - 

6FDmf 

1=NP 

S 

wP,Aco~Al 

n= co 

Zexp  (-n27r2fpUo) 

v  Pm yPc,o  J 

—2 

1=1 

a  . 

pi,o 

n= 1 

,[A] 


(36) 


In  Eq.  (36),  Np  =  total  number  of  size  ranges  of  the  active  material  particles  present  in  the 
cathode,  api  o  =  surface-area  fraction  mean  particle  radius  of  the  active  material  particles  of  the 

Dt 

particle  size  range  i ,  and  r  .  o  =  — —  (dimensionless  time  associated  with  a  .  o ). 

a 

pi,o 


The  total  amount  of  charge  stored  in  the  cathode  active  material  particles  in  time  t  during  the 
cell  discharge  period  is  given  by 


Q’  i 

Q‘  -Q'0  =  \dQt  =\V'cdt 

Q'o  0 


(37) 


The  substitution  for  IuC  from  Eq.  (36)  into  Eq.  (37)  and  the  integration  of  the  resultant 
expression  leads,  on  simplification,  to: 


A<2'=0'-e:  = 


f  ,r  3 

Fm'f 

i=Np 

Z 

v  Pm yPc,o  y 

i= 1 

wPiAco-c iX  1 


i-4z 

n= 1 


f  -n27T2Tnlr  ^ 


J) 


[coulomb]  (38) 


ll 


For  the  situation  of  (ca  -q )  being  independent  of  the  active  material  particle  size  distribution, 
i.e.,  (co-c,).  =  (c0-q),  Eqs.  (36)  and  (38)  lead  to: 


jt,C  _ 


r  6  FDm'f'' 


Pm 


r  iviyPc,0  J 


(Co~F) 


i=Nn 


r ...  \ 


w 


pi,o 


—2 

V  api>°  J 


n= oo 

^  exp  (~n27r2Tpio) 


,[A] 


(39) 


a<2' =G'-e:  = 


Em''c(co-Ci) 


Pi 


M  Pc, o 


i=N„ 


W  ■  < 

pi,o  ' 


^  «=1 


6  —rp'Tp'T  .  M 


,  [com /cm/?]  (40) 


Also,  for  the  active  material  particles  of  only  one  size  present  in  the  cathode,  Eqs.  (39)  and  (40) 
lead  to: 


I,c  = 


l.C 


6FDmo 

y  pMyPc.o  j 


c„  -c. 


V  ap-°  J 


^exp(-nVrpo),[A] 

n= 1 


(41) 


A<2r=e'-e:  = 


Frrif(co-cx) 


P M  yPc,o 


:=co 

\ 


i-A> 


V 


n 


,  [cow/om^]  (42) 


where  r  =  —  . 

a 

p,o 

It  was  experimentally  shown  [2]  that  the  mass  diffusivity  of  lithium  in  a  metal  phthalocyanine 
(e.g.  FePc)  cathode  is  a  function  of  both  the  stored  or  intercalated  lithium  concentration  and 
temperature.  Therefore,  the  diffusivity  of  lithium  in  a  MyPc  cathode  particle  is  given  by 

D(x,T)  =  Do(x)e~E°Dl(RT]  (43) 


where  Do  (x)  =  pre-exponential  factor  function  of  the  stored  lithium  content, 

f  .  .  A 


x  = 


moles  of  the  stored  lithium 


g-mole  of  the  MyPc  in  the  absence  of  stored  lithium^ 


Ea  D  =  activation  energy  for 


lithium  diffusion  in  MyPc  cathode  active  material  particles;  R  -  universal  gas  content - 
8.314 J  moFxK~x  and  T  -  absolute  temperature  . 


In  Ref.  [2]  it  has  been  shown  that  D(x,T)  at  a  temperature  (for  example,  20°C)  decreases  as  x 

increases  for  lithium  storage  in  the  active  material  FePc  particles.  The  mass  diffusivity  of  lithium 
averaged  over  the  stored-lithium  in  the  cathode  active  material  at  any  temperature  is  given  by  the 
following  expression 
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-E*.d/{RT) 


(44) 


D  = 


J  D0  (x)dx 

o _ 

X 


e 


f 

D  =  Do exp 

v 


-E  ^ 

‘-'a.D 

RT 


(45) 


Under  the  assumption  of  the  cell  operation  controlled  by  the  solid-state  diffusion  of  lithium  in 
the  cathode  active  material  of  MyPc,  the  voltage  loss  associated  with  the  lithium  ion  transport  in 
the  electrolyte  and  that  associated  with  the  electron  transport  in  the  acetylene  black  carbon  may 
be  assumed  negligible.  Simple  estimates  of  these  voltage  losses,  however,  can  be  obtained  from 
the  following  mathematical  expressions  at  any  time  during  the  cell  discharge  period.  Voltage  loss 
associated  with  the  ion  transport  in  the  electrolyte  present  in  the  composite  cathode  of  thickness, 
l  ,  can  be  estimated  from: 


AO 


c 

electrolyte 


CiC 


i:l 


r\  C— electrolyte 
AKeff-ion 


(46) 


where  the  effective  ionic  conductivity  of  the  electrolyte  in  the  cell  cathode,  V,/  el"i'lohu‘ ,  is  given 
by 


C -electrolyte 
K-eff-ion 


C- electrolyte 
Kion  £ 


C— electrolyte 


T 


C— electrolyte 


(47) 


In  Eqs.  (46)  and  (47),  /,(  =  cathode  electrode  geometric  current  density 


lc  =  thickness  of 


the  cathode  electrode,  KCum  e,ecrro,yte  -  ionic  conductivity  of  the  electrolyte,  sc  elec,roly,e  =  volume 
fraction  of  electrolyte  in  the  cathode,  and  Tc~elec,roly,e  -  tortuosity  factor  of  the  lithium 
ion-transport  pathways  in  the  electrolyte,  It  C  =  total  current  of  the  cell  cathode,  and  A(;  = 

cathode  geometric  area  in  contact  with  the  cell  electrolyte  separator.  The  factor  2  in  the 
denominator  of  Eq.  (16a)  appears  because  the  ionic  current  in  the  electrolyte,  at  any  time,  at  the 
cathode  electrolyte-cathode  current  collector  interface  is  zero  and  that  at  the  cell 
separator-cathode  interface  is  .  The  voltage  loss  associated  with  the  electron  transport  in  the 
carbon  black  present  in  cathode  can  be  estimated  from 


AO 


C-carbon 

electron 


CiC 


i:i 


?K-C 

^  ejf  -electronic 


(48) 
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where  ^fr -electronic =  effective  electronic  conductivity  of  cathode  carbon  and  is  given  by 


K 


^C— carbon 


eff  —electronic 


carbon  C— carbon 

v electronic _ 

^C— carbon 


(49) 


where  =  electronic  conductivity  of  the  carbon  black  present  in  the  cell  composite 

cathode,  sc~carbon  —  volume  fraction  of  the  carbon  black  in  the  cathode,  and  Tc-carbon  —  tortuosity 
factor  of  the  electron-transport  pathways  in  the  carbon  black.  It  is  here  noted  that  the  ratio  of  the 
LioPc  electronic  conductivity  [2]  to  that  of  the  carbon  black  is  of  the  order  of  10  10 — 10-11; 
indicating  that  the  electronic  transport  in  the  cell  cathode  takes  place  mainly  via  the  carbon  black 
present  in  the  cathode. 

4. 1.1.1  Reversible  or  Open- Circuit  Potential  of  the  Electrochemical  Reaction  at  the  Cell 
Cathode: 


The  electrochemical  reaction  that  takes  place  at  the  boundary  of  three  phases;  a  solid  electronic 
conductor  (e.g.  carbon  black),  cathode  active  material  (e.g.  LixMyPc),  and  an  electrolyte  (e.g. 
LiBF4  in  a  microporous  poly(ethylene  oxide));  is  represented  as 


Li 


(in  electrolyte) 


+  e‘ 


(in  carbon  black) 


LixMyPc 


x  g-mole  of  the  intercalated  lithium 
present  in  1  g-mole  of  initial  MyPc 


□  Lix+1MyPc 


(x+l)  g-mole  of  the  intercalated  lithium 
present  in  1  g-mole  of  initial  MyPc 


(50) 


The  intercalated  lithium  is  treated  as  a  solute  present  in  solid-state  solvent,  e.g., 
MyPc  [5,  6,  7]. 

The  thermodynamic  equilibrium  condition  for  the  reaction  (50)  is 


echem 


echem 


JU  .  .  +  JU 

r  Li  (in  solid  electrolyte)  '  e  (in  solid  carbon  black) 

chem  .  chem 

XMlA  (in  solid  MyPc)  ~l~  MvLPc(solid) 


=  (*  +  !  )ju£em 


chem 


(in  solid  MyPc)  /4lyPc(solid) 


(51) 


Equation  (51)  leads  to 


echem  _l_  echem 

'  Li+  (in  solid  electrolyte)  '  e"  (in  solid  carbon  black) 


_  chem 

—  Mi 


(in  solid  M  Pc) 


(52) 


Expressing  the  species  electrochemical  potentials,  jU(‘ch,'m  in  Eq.  (52),  in  terms  of  chemical  and 
electrostatic  potentials,  leads  to 

chem  .  /^(T) 

_ '  Li+  (in  electrolyte)  +  ^Li+  ^  electrolyte 


+ 


chem  , 

H  ,  .  +  z  - 

1  e  (in  carbon)  e 


F  O 


carbon 


chem 

’  "u  (in  M  Pc) 


(53) 


Simplification  leads  to: 
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f(® 


—  (1) 

electrolyte  carbon 


)- 


chem  chem  chem 

/^Li(in  My  Pc)  '  Li+  (in  electrolyte)  '  e_  (in  carbon) 


(54) 


Equation  (54)  is  rewritten  as 


(® 


carbon  ^  electrolyte 


)- 


chem  .  chem  chem 

*  Li+  (in  electrolyte)  '  e~  (in  carbon)  '  Li(in  MyPc) 


(55) 


where  ® 


electrolyte 


—  electrostatic  potential  in  the  electrolyte  just  outside  of  the  double  charge  layer 


carbon  ' 


at  the  three  phase,  carbon  black-cathode  active  material-electrolyte  boundary  and  ® 
electrostatic  potential  in  the  carbon  black  electronic  conductor  in  the  cell  cathode.  By  definition, 
( c[/  -  d>  ,,,,)  =  Ec  =reversible  potential  of  the  electrochemical  reaction  that  occurs  at  the 

\  carbon  electrolyte  I  rev  a 

cell  cathode  for  the  cell  open-circuit  conditions.  Now,  Eq.  (55)  is  written  as 


Ec  = 

rev 


chem  chem 

‘  Li+  (in  electrolyte)  '  e~  (in  carbon) 


chem 

"  /^Li(in  M  Pc) 


(56) 


chem 


Expressing  ^ electrolyte}  and  //u(inMyPc) 


as 


chem  o  ,  n'T'  1  „ 

juT.+/.  ,  ,  =jUT  +  ,  ,  +Ri  lnaI  +  , 

‘  Li  (in  electrolyte)  ‘  Li  -electrolyte  Li  -electrolyte 


(57) 


and 


EnZrc 


Li+-M„Pc  +  RT  1,1  aLi-MvPc 


(58) 


where  ^  (  =  standard-state  chemical  potential  of  Li+  in  the  solid  electrolyte,  e.g.,  when 

CLi+ -electrolyte  =  ^  4  A7 ;  //"  M  Pc  =  standard- state  chemical  potential  of  lithium  in  the  cathode  active 
material,  e.g.,  when  the  intercalated  lithium  concentration,  c°u_M  Pc  =  x°cM  Pc  0 ,  where 

x"  =  10 ~9  g -mole  of  the  intercalated  lithium  per  g-mole  of  MyPc  and  cM  Pc  o  =  initial  molar 

concentration  of  the  cathode  active  material  MyPc;  the  activity  of  Li+  in  the  electrolyte, 

c  .+ 

a  +  =  V  Li  -electrolyte  an(j  acqvjty  Qf  lithium  in  MyPc, 

Li  -electrolyte  '  Li  -electrolyte  J  j 

Li+  -electrolyte 


a 


Li-M„Pc 


=  Yu- 


CLi-M„Pc 


MPc  o 


CLi-M„Pc 


■  =  Yu- 


-M„Pc 


xc\ 


X  c„ 


=  Tu¬ 


rn  „Pc 


yx°J 
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Inserting  the  information  from  Eqs.  (57)  and  (58)  into  Eq.  (56)  leads  to 


El  = 


U  .  +  jU  _ 

'  Li  -electrolyte  '  e  -carbon 


"  /Li-MyPc 


RT 


In 


a  , 

Li  -electrolyte 
V  flLi-MyPc  j 


(59) 


Note  that  in  obtaining  of  Eq.  (59),  the  actual  state  of  a  charge  conducting  electron  in  carbon  has 
been  taken  to  be  equal  to  its  standard  state.  Equation  (59)  is  rewritten  as 


Ec  =EC’°  + 

rev 


f  RT\ 


V  r  J 


In 


y 

'  Li+ -electrolyte 


Li+  -electrolyte 

.  c°.+  . 

I  Li  -electrolyte  J 


Yu- 


M.,Pc 


v*  ) 


(60) 


where  the  standard-state  electric  potential  of  the  reaction  that  occurs  at  the  cathode,  Ec'°  is 
given  as 


Ec'°  = 


^U+ -electrolyte  -carbon  Ru  MyPc 


(61) 


Equation  (60)  is  re-written  as 


7 c  =  EC-° 

^ rev 


f  RT ' 

ln< 

+ 

V 

r, 


Li+  -electrolyte 


A 


Li+  -electrolyte 


y  Li+  -electrolyte  J 


f  RJ\ 


V  r  J 


ln  \Yu- 


\X  J 


(62) 


Equation  (62)  can  be  written  as: 


Ec  =El'°  + 

rev 


C,o 


r  RT\ 


V  r  J 


In 


y  c 

'  Li+  -electrolyte  Li+  -electrolyte 


Li+  -electrolyte  J 


-In 


^  T'n-MyPc  ^ 
V  J 


f  RT  ^ 


\  r  J 


lnx  (63) 


Ec  =Er 

rev  G 


C,o 


f  RT  \ 


V  r  J 


lnx 


(64) 


where 


EC,°  =  EC,o  + 


RT 

\~y  j 


In 


y  c 

'  Li+ -electrolyte  Li+ -electrolyte  , 

*  ~  *  JC 

(  ^Li-MyPc  Li+ -electrolyte  J 


(65) 
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For  a  given  electrolyte  at  a  temperature,  E((y"  would  be  constant  with  respect  to  x  if  variation 
in  yu  M  Pc  is  assumed  to  be  negligibly  small.  Equation  (64)  suggests  that  Ecnv  has  a  linear 
relation  with  In  x  •  provided  yu _M  Pc  is  constant  with  respect  to  x.  So,  the  plot  of  measured 


values  of  Ef  vs.  lnx  should  turn  out  to  be  a  straight  line  with  slope  = 


-RT 


.  If  so,  then, 


v  F  j 

one  can  determine  ECG°  from  the  measured  values  of  Ecm  vs.  x  using  Eq.  (64).  Finally,  the 
relation,  Eq.  (64),  can  be  employed  to  compute  Ecrev  as  a  function  of  x  for  the  prediction  of 
performance  of  a  lithium-ion  cell. 


Equation  (62)  can  also  be  expressed  as 


Ec  =EC’°  + 

rev 


' RT ' 

r 


In 


y 

'  Li+ -electrolyte 
v  y Li-MyPc  j 


f  RT\ 


V  r  j 


In 


c  * 

Li  -electrolyte 
C 

\  Li + -electrolyte  J 


f RT V  (  x^ 

In 


V  F  J 


\X  J 


(66) 


For  the  assumption  of 


v 

'  Li + —electrolyte 
v  y Li-MyPc  j 


=  constant,  Eq.  (66)  suggests  that  Efev  decreases  as  x 


increases  above  x°  for  a  fixed  c  +  ,  ,  ;  whereas  increases  as  c  ,+  ,  ,  increases 

Li  -electrolyte  rev  Li  -electrolyte 

above  c°  .  for  a  fixed  x  value  at  a  temperature. 

Li  -electrolyte  A 


4. 1.1. 2  Electrode  Kinetics  at  the  Cell  Cathode: 

The  cathode  electrochemical  reaction  takes  place  at  the  boundary  of  three  phases  -  a  solid 
electronic  conductor  (carbon  black),  cathode  active  material  and  an  electrolyte.  The 
electrochemical  reaction  is  represented  as 


Li  ,  ,  +  c 

(in  electrolyte) 


(in  carbon  black)  hi  y  PC  ^ 


{ in  Ig-mole  of  the  solid-'] 
state  solvent  of  MvPc 


(67) 


It  is  here  assumed  that  the  cathode  active  material,  MyPc,  is  a  solid-state  solvent  for 
lithium  storage  of  the  intercalated  lithium  during  the  period  of  cell  discharge.  The 
forward  or  cathodic  direction  process  of  the  reaction,  Eq.  (67),  which  has  units  of 


g  -  mols  of  the  forward  reaction  process 
time  □  active  material  interf acial  area 


,  is  simulated  by 


r;  =k 


c,  1 


^I.i 1  in  electrolyte^MyP(T 


-pcEcFj(RT) 


(68) 


where  J3C  =  symmetry  factor;  indicative  of  the  fraction  of  the  voltage  difference,  Ec  , 
(difference  between  the  electrostatic  potentials  at  the  three-phase  boundary  and  in  the  electrolyte 
just  outside  the  double  charge  layer)  involved  in  the  cathodic  direction  charge  transfer  process; 
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a  .  —  activity  of  lithium  ions  in  the  electrolyte  = 

Li  -in  electrolyte  J  J 


r, 


Li+  ,conc 


c° 

V  Li+  J 


(dimensionless), 


/  +  nc  =  activity  coefficient  of  lithium  ions  in  the  polymer  electrolyte,  c[j+  =  molar 

concentration  of  lithium  ions  in  the  electrolyte,  c"j+  =  molar  concentration  of  Li+  in  the  polymer 
electrolyte  when  the  electrolyte  active  component,  e.g.,  LiBF4  is  present  in  ‘infinitely’  dilute 
form,  e.g.,  c°  +  =  10  4 M  ,  taken  to  be  the  standard  state  of  Li+  in  the  electrolyte  at  which  state 
Y  ,+  =1  and  a  ,+  .  ,  =  I ;  a,, ,,  =  activity  of  the  cathode  active  material,  MvPc,  here 

'  Li  ,conc  Li+ -in  electrolyte  MyFc  J  1 


a 


M.,Pc 


may  be  taken  equal  to  1  because  MyPc  is  assumed  to  be  a  solid  solvent  for  lithium 

storage;  F  -  Faraday’s  constant  =  96487  coulomb  of  charge  per  g-equivalent  or  g-mole  of 
‘unitly’  charged  particles  (e.g.,  electron,  lithium  ion);  R  =  gas  constant  = 

8.314  J  mol'1  K1  ;T  =  temperature  at  the  three-phase  boundary,  [AT];  and  kcf']  =  rate  coefficient 


for  the  forward  or  cathodic  direction  process  of  the  reaction,  Eq.  (67), 
equation  (68)  is  rewritten  as 


g-mole 


cm  -sec 


.  The  rate 


ri=\Kra^Ju-conc 


) 


VCLi+  J 


~/3c  ECF /(RT) 


(69) 


=  k. 


CLi+ 

c° 

V  Li+  J 


-Pcecf/(rt) 


where 


(70) 


4  =(4J 


VI y Pc j  _conc 


g-mol 
cm 2  -  sec 


(71) 


The  rate  of  the  anodic,  oxidation  direction  process  of  the  reaction,  Eq.  (67),  is  assumed  to  be 
dependent  on  the  intercalated  lithium  activity  in  the  cathode  active  material  MyPc.  The  anodic 
direction  reaction  process,  with  units  of 

g -moles  of  the  anodic  direction  reaction  process  occuring 


cm  ■  sec 


,  is  represented  by 


—  JrC1 
~  Kb  a,Li-inM„Pce 


(72) 


where  (l  -  j3c  )  =  fraction  of  the  voltage  difference,  Ec  ,  involved  in  the  anodic  direction  charge 
transfer  process;  oLi.inM  Pc  =  activity  of  lithium  intercalated  in  the  active  material  MyPc 
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=  Yu, 


XC, 


cone  o 

CLi-in  MvPc 


=  Yu, 


cone  o 
X  C 


=  Yu, 


f  x  ^ 


Yu  ,  cone 


=  activity  coefficient  of  lithium  in  the 


MvPc 


V*  J 

=  molar  concentration  of  the  intercalated  lithium  in 


cathode  active  material,  MyPc,  cLi_inM  Pc 
MyPc,  c°i  in  M  Pc  =  molar  concentration  of  lithium  in  MyPc  when  it  is  present  in  MyPc  in  the 
‘infinitely’ dilute  form,  e.g.,  x"  =  10  9,  for  which  state  «Li.inMyPc  =1,  7LiiConc  =1,  x=  g-mole  of 

the  intercalated  lithium  per  g-mole  of  MyPc;  k(h  ']  =  reaction  rate  coefficient  for  the  anodic 
direction  charge  transfer  process.  Equation  (72)  is  rewritten  as 


rbC  =(tf’lYu,  cone) 


V*  J 


(1 

RT 


(73) 


t{  =  ki 


f  x^ 


\X  ) 


(1-/T  j  EXF 
RT 


where, 


(74) 


g =(gJrLm). 


g-mole 


cm  -sec 


(75) 


is  the  rate  coefficient  for  the  backward  or  anodic  direction  charge  transfer  process. 

The  net  reaction  rate  at  the  cathode  during  the  period  of  lithium  intercalation  into  the  cathode 
active  material,  i.e.,  during  the  cell  discharge  period,  is  given  by 


f 

r 

y 

-PCECF 

fx 

\ 

(l  -pc)ecf 

\rc\-  rc  -rc  - 

kc 

Li+ 

p  RT 

-kC\ 

\p  RT 

r  \rf  r>’  ~ 

Kf 

Kb  \ 

o 

lCLi+ 

) 

\X 

) 

(76) 


Because  c°.+  and  x°  are  constants  for  the  selected  standard  states,  Eq.  (76)  is  transformed  to: 


rc  = 


k„c  +e 

ff  Li+ 


-P  E'-F 
RT 


kbbxe 


i -pc')EcF 
RT 


(77) 


where  kcff 


(  kc  > 
Kr 

cm 

and  kchb  = 

(kc 

Kb 

g-mol 

o 

5 

o 

2 

lCLV  J 

_sec^ 

l X  J 

cm  -sec 

The  net  reaction  rate,  | | ,  is  related  to  the  interfacial  charge  transfer  current  density,  |zsc  | ,  at  the 
cathode  by  the  following  equation: 
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(78) 


-pL  elf 


(l -fic')EcF 


k„c  ,e 

ff  Li+ 


RT 


~hbxe 


RT 


Or, 


L-C 


=  r 


s\_kc 

~Kff 


(1  -Pc)ecf 


xe 


CLi+e 


Kc 


(79) 


kc„ 


where  Kc  =  — ,  [, g-mole  cm  3  is  the  kinetics  based  equilibrium  constant  for  the  reaction  at 


xbb 


the  cell  cathode. 


For  the  reaction  equilibrium  at  the  cathode,  =  0 ,  where 


Ec  =EC 

rev 


(80) 


Using  Eq.  (79)  in  Eq.  (80)  and  simplification  leads  to 


preq  p  RT  _ 
CLi+ 


Fxqe 


J  E,evF 
e<l„  RT 


Kc 


(81) 


where  xeq  is  the  number  of  g-moles  of  the  intercalated  lithium  per  g-mole  of  the  cathode  active 
material,  MyPc,  for  the  equilibrium  between  the  electrolyte  and  the  active  material  with  respect 
to  lithium. 

Equation  (81)  is  re-expressed  as 


Vc req  p  RT 


=  Fkbbxeq' 


RT  — 


=  l 


;C 


(82) 


where  f  =  exchange  current  density,  [A-cmJerface  J  and  Ecrev  -  potential  difference  between 

the  electrostatic  potential  in  the  electronic  conductor  at  the  three-phase  boundary  and  that  in  the 
electrolyte  just  outside  the  double-charge  layer  at  the  equilibrium  state. 

Determination  of  Ec  and  Ecrev  is,  in  general,  difficult.  So,  it  is  here  suggested  that  Ec  and 
Ecm  may  be  measured  relative  to  a  reference  electrode.  The  surface  overpotential  at  the  three 
phase  boundary  of  the  cell  cathode,  due  to  a  finite  value  of  I  if  I ,  is  defined  as 
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(83) 


rf  =EC -Ec 

•  s  rev 


The  surface  overpotential,  i][  ,  is  also  called  the  activation  polarization  voltage  change 
associated  with  the  electrochemical  reaction  involving  the  charge  transfer  process  across  the 
cathode  three  phase  boundary.  Substituting  Ec  =  Ecrev  +  rf  in  Eq.  (78)  leads  to 


ic  = 


-P  Ere,F 

kcffFc  +e  RT  t 

ff  Li+ 


KbF  xe 


> -Pc)eLf 

RT 


(l -0C)r,fF 


(84) 


C 
s 

The  cathodic  charge  transfer  coefficient,  acc  ,  and  the  anodic  charge  transfer  coefficient,  aca  , 

are  defined  as  follows :  of  =  /3cn  =  f3c ;  acu  =  (l  -  J3C  )  n  =  1  -  J3C ;  because  for  the 

electrochemical  reaction  at  the  cell  cathode,  Eq.  (67),  n,  the  number  of  electrons  transferred  per 
one  reaction-event  is  equal  to  1. 

The  expression  for  the  current  density  at  the  cathode  three  phase  boundary  from  Eq.  (85)  is  now 
written  as 


i 


=  f 


\  -PCrfsF 


Li 

req 

VCLi+  J 


f  X  ^ 

\xeq  p 


-Pc)nfF 


(85) 


f  \ 

r 

C  C  r. 

-ac  rjs  F 

f  X  >1 

C  Cr 

a  a  rjs  F 

II 

CLi+ 

P  RT  _ 

\p  RT 

VLi+  J 

(86) 


Equation  (86)  suggests  that  as  x  increases  during  the  cell  discharge  period,  the  current  density 
should  decrease.  When  the  intercalated  lithium  in  the  cathode  material  attains  a  value  of  j  =  xeq , 
rf'  =0,  and  cT.+  =c*f+;  |ic|  becomes  equal  to  zero. 

* A  Li  Li  |  s  |  x 


Equation  (82)  is  used  to  obtain  the  following  equation 


kc  req 
KffCu* 

kC  Yeq 
F-bbX 


(87) 


Taking  log?  of  both  sides  of  Eq.  (87)  leads  to 


Ec  = 

rev 


F  J 


Hb  J 


(88) 
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'  RT' 

f  kc'\ 

f  RT 

( req  ^ 
*-'T  •  + 

in 

JJ 

+ 

in 

Li 

,~F  j 

kc 

yKbb) 

Cf  j 

Kx  ) 

Or,  Ec  = 

’  rev 


Equation  (88)  is  the  desired  formula  to  compute  Ecm  based  on  the  kinetics  of  the 

electrochemical  reaction  that  occurs  at  the  cell  cathode.  Note  that  Eq.  (88)  is  similar  to  Eq.  (64) 
derived  using  the  thermodynamic  principle  of  reaction  equilibrium  to  compute  the  reaction 
equilibrium  electric  potential  at  the  cathode  active  material -electrolyte  interface  for  a  given 
electrolyte  and  cathode  active  material,  such  as  MyPc  at  a  cathode  temperature. 


If  one  knows  the  experimental  values  of  Ecrev  and  xeq ,  one  can  determine 


fkc\ 

kc 


for  a  given 


electrolyte  from  Eq.  (88).  Inserting  the  information  given  in  Eq.  (87)  into  Eq.  (82)  leads  to 


(89) 


Notice  the  relation  of  the  exchange  current  density,  if ,  with  the  chemical  and  electrochemical 
kinetic  parameters  without  the  involvement  of  the  reversible  electrochemical  reaction  potential, 
E(m, ,  at  the  cell  cathode. 


From  Eq.  (86),  with  c  +  =  c"'  ,  particularly  for  the  situation  of  the  lithium  diffusion  in  the 
cathode  active  material  particles  controlling  the  overall  cell  performance,  one  obtains: 


( 

(af+acc)vff 

-ccqnfF 

II 

o 

1- 

X 

je  RT 

e  RT 

1  1  ° 

Y661  , 

yx  J 

(90) 


Taking  log?  of  both  sides  of  Eq.  (90)  and  further  simplification  leads  to 


In 


1- 


exp 


RT 


(91) 


One  should  note  that  ( r/;  +  c/  j 
(67). 


(l-  j3c^jn  +  j3cn  =  n  =  1  for  the  electrochemical  reaction,  Eq. 
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Equation  (91)  suggests  that  In 


1- 


f  x  3 


\X  J 


exp 


RT 


)€F 


plotted  versus  (-?7.sc)  should 


be  a  straight  line  with  slope  = 
(—rf)  is  varied. 


f  acF^ 


v  RT  j 


and  intercept  =  In  rf  ,  provided  rf  does  not  vary  as 


Such  a  plot  based  on  the  experimental  data  of  |/(;  | ,  x,  x" ,  and  (-rf  j  should  be  used  to  obtain 
electrokinetic  parameters  rf  (=  J3cn)  and  rf ,  the  exchange  current  density  for  the 
electrochemical  reaction,  Eq.  (67). 

For  the  low  surface  overpotential,  (—rf  )  ;  i.e.,  for  a  low  value  of  |/f  | ,  the  Taylor  series 

-gfi/fF  acarfsF 

expansion  of  e  RT  and  e  RT  in  Eq.  (86)  leads  to 


rf  =  rf 


f  r  S  f 

1-  — 


V 


xeq  j 


C  C 

a ;  +a, 


v 


Vxe“jj 


(92) 


For  the  cathode  active  material  composition  state  close  to  that  in  equilibrium  with  the  electrolyte, 
i.e.,  for  x  =  xeq ,  Eq.  (92)  leads  to 


i;  |  =  i 


ic  (ac 

o  y  c 


(93) 


This  equation  suggests  that  the  experimental  data  on  |/(|  obtained  as  a  function  of  (—rf)  for 
low  values  of  (-rf) ,  with  the  cathode  active  material  composition  state  close  to  its  state  in 
equilibrium  with  the  electrolyte,  can  be  employed  to  obtain  rf  experimentally. 

For^af  (-rf  )  F)  »  RT ,  Eq.  (86)  reduces  to 


Or, 


■c  ■ c 
i,  =ine 


-a£  r/s  F 
RT 
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RT 

rfF 


\  fi.-ch 

In 
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■c 
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l 


J 


(94) 


(95) 
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For  acarfsF  »  RT ,  Eq.  (86)  reduces  to 


ic  =  -if 


f  x  \ 


\xeq  J 


RT 


(96) 


Or, 


In 


=  In 


f  v  \ 


\xqj 


v  RT  j 


Vs 


(97) 


Equation  (97)  leads  to 


Vs  = 


r  RT  A 

C 


In 


-i 


c  A 


•c 


V  *  y 


(98) 


Plot  of  ?/5  vs.  In 
slope  is  either 


fi^ 


J 


or  plot  of  r/s  vs.  In 


H, 


x 


eq 


o  \  X  j 


r  -RT^ 


\acF  j 


or 


r  RT  ^ 


\  aaF  J 


is  called  the  Tafel  plot.  The  Tafel 


4.1.2  Electrolyte  Separator  Formulation 

The  geometric  current  flow  density  through  the  electrolyte  separator  of  microporous 
poly(ethylene  oxide)  containing  LiBF4  is  computed  from 


,t,c 


if  =  ■ 


(99) 


where  Af  =  geometric  area  of  the  separator  perpendicular  to  the  current  flow  direction.  For  the 

situation  of  the  lithium-ion  cell  performance  controlled  by  the  diffusion  of  intercalating  lithium 
in  the  cathode  active  material,  the  lithium  ion  concentration  gradient  in  the  polymer  separator  is 
assumed  to  be  negligible.  The  voltage  loss  in  the  electrolyte  separator  is  given  by 


AO 


s 

electrolyte 


;S7S 


i.l 


S  -electrolyte 
Keff 


(100) 
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where  Is  =  thickness  of  the  electrolyte  separator  and  KSeff elecrrolyte  =  ionic  conductivity  of  UBF4  in 
microporous  poly(ethylene  oxide). 

4.1.3  Anode  Formulation 

4.1.3.1  Electrode  Kinetics  at  the  Cell  Anode 

For  the  anode  of  a  thin  lithium  sheet  in  contact  with  the  electrolyte  separator,  the  electrochemical 
reaction  occurring  at  the  anode-electrolyte  interface  is  represented  as 

kA1 

Li(s)<===±Li+  (in  electrolyte)  +  c  (in  metal  of  Li)  (101) 

K ' 


The  forward  reaction  rate  of  the  electrochemical  reaction,  Eq.  (101),  which  has  units  of 
g-mole 


sec- cm'  interface 


is  expressed  as 


A  1  A,1 
rf  =kf  au (sf 


i-paY 


=ky-n‘ 


(102) 


where  kAA  =  forward  reaction  rate  coefficient, 


g-mole 

sec-  cm2  interface 


a|  i(  v)  =  activity  of 


solid-state  lithium=l  (a  usual  assumption),  /3A  =  symmetry  factor  for  the  electrochemical 

reaction,  Eq.  (22a),  EA  =  difference  between  the  electrostatic  potential  in  metal  of  lithium  and 
that  in  the  electrolyte  just  outside  the  double-charge  layer  at  the  cell  anode-electrolyte  interface, 
[V ] ,  when,  for  example,  the  cell  is  being  discharged. 


The  backward  direction  reaction  rate  of  the  reaction,  Eq.  (101),  which  has  units  of 
g-mole 


sec  -cm2  interface 


,  is  represented  as: 


A  ,  A.l  -BaEaf/rT 

r.  =  k.  a  ,+/.  ,  ,  je  H  1 

b  b  Li  (in  electrolyte) 


(103) 


where  kAA  =the  backward  direction  reaction-rate  coefficient  for  the  reversible  reaction,  Eq. 


(101), 


g -moles 


sec- cm'  interface 


a . 


cu 


=  y  ,+  =  activity  of  Li  in  the  electrolyte  solution 


CLC 


adjacent  to  the  anode  electrode-electrolyte  interface  just  outside  the  double-charge  layer,  cu+  = 
molar  concentration  of  Li+  in  the  electrolyte  solution  adjacent  to  the  electrode-electrolyte 
interface,  and  c".  =  standard- state  concentration  of  Li+  when  it  is  present  in  ‘infinitely  dilute 

form’  in  the  solvent,  here,  poly( ethylene  oxide),  yu+  =  activity  coefficient  of  lithium  ions  in  the 
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electrolyte  to  account  for  the  non-ideal  behavior  of  Li+  in  the  solvent.  For  Li+  ions  present  in  an 
‘infinity  dilute  form’  in  the  electrolyte,  e.g.,  c"  =  10  4A7,  /  =  /"+  =  1  and  a'.  =  1 . 

Equation  (103)  is  rewritten  as 


r?  = 


Kr^) 


c° 

V  CLi+  J 


-Paeaf/rt 


rA  —  kA 

'b  —  Kb 


c° 

V  Li+  J 


-Paeaf/rt 


where 


(104) 


(105) 


(106) 


and  has  units  of 


g-mole 

sec- cm2  interface 


The  reaction  rate  expressions,  Eqs.  (102)  and  (104),  are  now  represented  as 


where 


A  (\-Pa)eafIrt 

V 


rA-kAr  ~Paeaf/rt 

'b  ~  KbbCu+t 


g-mole 

sec-  cm2  interface 


( kAAv 

Kb  ^Li+ 

cm 

c° 

V  Li+  J 

_sec_ 

(107) 

(108) 


(109) 

(HO) 


The  net  reaction  rate  at  the  anode  during  the  cell  discharge  period,  which  has  units  of 
g -moles  of  reaction  (101)  occuring 


sec-  cm  interface 


is  given  by 


_  La  I  _  A  _  A  _  rA  (i-Pa)eaf/rt  _  A  -paeaf/rt 
F  -y  \~ rf  rb  - Ktfe  Kbbe  cLi+ 


(111) 
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where 


iA  |  =  interfacial  charge-transfer  current  density, 
discharge. 

For  the  reaction,  Eq.  (101),  to  be  at  equilibrium; 

|C|  =  0  (112) 


cm 


during  the  period  of  cell 


Using  the  result  from  Eq.  (1 1 1)  in  Eq.  (1 12)  leads  to 


Fkffd'1~^E™F/RT  =  Fk*be~fiAE™F/RTcZ  =  iA 


(113) 


where  iA  =  exchange  current  density  for  the  electrochemical  reaction,  Eq.  (101);  and  E*ev  = 

difference  between  the  electrostatic  potential  in  the  anode  electrode  and  the  electrostatic  potential 
in  the  electrolyte  adjacent  to  the  electrode-electrolyte  interface  for  the  anodic  reaction,  Eq.  (101), 
to  be  at  equilibrium. 

The  surface  overpotential,  with  units  of  volts,  at  the  electrode-electrolyte  interface,  rjf ,  is 
defined  as 


(114) 


The  surface  overpotential,  r/f  ,  is  also  called  the  activation  polarization  voltage  change 
associated  with  the  electrochemical  reaction,  Eq.  (101),  involving  the  charge  transfer  process 
across  the  electrolyte-anode  electrode  interface.  Inserting  EA  =  ( EAv  +  rjf  j  into  Eq.  (Ill)  leads 

to 


if  = 


kffFe 


'l 

RT 


(l -pA)Fr,f 
RT 


-kbbFcu*e 


~PAFE%V  -p  Frjf 
RT  RT 


(115) 


Using  the  information  from  Eq.  (113)  in  Eq.  (115)  leads  to 


(l -pA)ntF 

f 

r 

-pAr,?F 

A 

II 

e  RT  - 

CLi+ 

e  RT 

r^q 
\  U+ 

J 

_Cmirf_ 

(116) 


Note  that  r/A  -  EA  -  EAv  is  a  positive  number  quantity  during  the  period  the  cell  is  providing 
electric  power  to  an  external  circuit. 

From  Eq.  (1 13),  one  finds  that 
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ErevF 

,  RT  — 
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kA 

V  Kff  J 


Take  loge  of  both  sides  of  Eq.  (117)  to  obtain 
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Ea  =  Ea'°  + 

rev  kinet 


f  RT\ 
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ru+cu+ 


V  Li+  J 


where  KL  = 


f  rt' 

r  K’1 ) 

in 

CT  j 

K 

V  ff  7 

(117) 


(118) 


(119) 


(120) 


Equation  (119)  is  the  kinetics  based  expression  to  compute  the  anode  reaction  reversible  electric 
potential  or  the  cell  open-circuit  anode  electrode  potential.  At  the  cell  open-circuit  condition, 


<?+  =  c  .  ,  lithium-ion  concentration  in  the  bulk  electrolyte.  Note  that  EAev  =  EA’°  = 


7  A,o 


standard-state  electric  potential  of  the  electrochemical  reaction,  Eq.  (101),  for  ceq+  =  c  +  =  c°+ . 


Equation  (119)  tells  us  that  as 
electric  potential  increases. 


(  eq  \ 

(  \ 

c  \ 

c  + 

Li 

Li+ 

c°  + 

c°  + 

V  Li  7 

^  Li+  7 

is  increased,  the  anode  electrode  reversible 


Using  the  information  from  Eq.  (117)  in  Eq.  (1 13)  leads  to 


(121) 
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Note  that  c‘J(  =  cu+  =  concentration  of  Li+  ion  in  the  bulk  electrolyte  solution.  The  anode 

electrode  reaction  exchange  current  density,  if ,  as  shown  by  Eq.  (121),  is  independent  of  the 
electric  potential. 

For  the  situation  of  the  intercalated  lithium  diffusion  in  the  cell  cathode  controlling  the  overall 
cell  performance,  it  is  likely,  c  +  =  cf.+ .  Equation  (116),  then,  reduces  to 


(l -PA)vfF 

-PAritF 

\iA\  =  iA 

I  s  |  o 

e  RT 

- e  RT 

(122) 


For  low  r/f  values,  i.e.,  for  low  if  values,  Eq.  (122)  is  reduced  to 


=  i 


RT 


Vs 


(123) 


This  equation  suggests  that  the  experimental  data  on  |/'4  as  a  function  of  r/f  can  be  employed 

to  experimentally  determine  the  anode  electrode  exchange  current  density,  if ;  for  its  application 
in  Eq.  (122).  Equation  (122)  can  also  be  expressed  as 


aAn?F  -aAnAF 

II 

€  RT  _e  RT 

(124) 


where  af  =  (4  -  J3A  )  n  =  (4  -  (3A  )  =  anodic  direction  charge  transfer  coefficient  across  the  cell 

anode-electrolyte  interface,  a f  =  (5An  =  fiA  =  cathodic  direction  charge  transfer  coefficient, 
(here,  n=l). 


For  large  current  density,  if  ,  values,  i.e.,  for  ( afrjfF )  »  (RT ) ;  Eq.  (124)  reduces  to 


k  I -aIN 


.•a 


a*n?F 

=  e  RT 


(125) 


Take  loge  of  both  sides  of  Eq.  (125)  to  obtain 


In 


'■  1 

l  *T  ) 

Vs 


(126) 
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This  equation  tells  us  that  a  plot  of  experimental  data  on  In 


versus  rjs  should  be  a 


straight  line  with  its  slope  = 


RT  ' 


This  could  help  one  to  determine  a A  experimentally. 


Because,  here,  for  one  electron  reaction,  (aA  +  aA  )  =  n  =  1 ;  once  aAa  is  known,  one  can 

determine  aA  .  These  experimentally  determined  charge  transfer  coefficients  aA  and  aA 
then  be  used  to  compute  the  anode  electrode  charge  transfer  current  density,  or,  the  anode 
electrode  polarization  voltage,  77  f  ,  for  a  current  density,  I/'4  from  Eq.  (124). 


can 


For  the  intercalated  lithium  diffusion  in  the  cathode  active  material  controlling  the  cell  overall 
behavior  during  the  cell  discharge  period,  the  anode  interfacial  charge  transfer  current  density, 
l'4 1 ,  should  be  computed  from 


(127) 


where  AA  =  geometric  area  of  the  cell  anode  electrode  sheet  in  contact  with  the  polymer 
electrolyte  and  |/f  C|  =total  cathode  current  as  predicted  by  Eq.  (39)  or  (41). 

4.1.3.2  Reversible  or  Open-Circuit  Potential  of  the  Electrochemical  Reaction  at  the  Cell 
Anode 

The  electrochemical  reaction  occurring  at  the  lithium  anode-electrolyte  is 

^  (in  electrolyte)  ^  (in  lithium  metal)  (128) 

The  thermodynamic  equilibrium  condition  for  the  electrochemical  reaction,  Eq.  (128),  dictates 
that 


chem  echem  ,  echem 

VuU)  =  M^  +  Me- 


Uchem  =  (  uchem  +  7  f  (T)  )  +  (  uchem  +  7  /  ([} 

'Li  ( .v )  ( 'l.i  (in  electrolyte )  Li+  electrolyte)  ( '<■  (in  Li  metal)  e 


A 

Li(r) 


(129) 

(130) 


This  equation  leads  to 


I?  A  _  ( ^A  _  ^ 

^  rev  y  Li  (s)  electrolyte  ^ 


f  chem  .  chem  chem  \ 

'  Li+  (in  electrolyte)  *  e~ {m  Li)  'Li (.s) 


(131) 


where  EAv  =reversible  electric  potential  of  the  electrochemical  reaction,  Eq.  (128),  occurring  at 
the  anode  electrode;  =  electrostatic  potential  in  the  lithium  metal  electrode,  Oelectrolyte  = 

electrostatic  potential  in  the  electrolyte  adjacent  to  the  anode  electrode  just  outside  the  double 
charge  layer;  electrolyte)  ’  ’  an<^  Ru(")  are  the  chemical  potentials  of  Li+  ion  in  the 

electrolyte,  electron  in  lithium  metal,  and  of  lithium  metal  of  the  cell  anode  electrode, 
respectively.  The  species  chemical  potentials  are  now  expressed  as 

uha;  ,  ,  =jur+RTlna.+l.  ,  ,  (132) 

'  Li  (in  electrolyte)  ~  Li  Li  (in  electrolyte)  v  7 
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(133) 


chem 

Ru(s) 


'■  Ru(s) 


+  RT\na 


Li  W 


(134) 


where  //+  =  standard-state  chemical  potential  of  Li+;  the  standard-state  of  Li+  ion  is  here  taken 
to  be  its  state  when  present  in  ‘the  polymer  solvent’  in  an  infinitely  small  amount,  e.g., 
c'j  =  10  4  A/  concentration  of  Li+  (or  of  LiBF4  in  the  polymer  electrolyte); 

c  \ 


a. 


Li+  (in  electrolyte) 


r r 


a 


(in  Li) 


=  electron  activity  in  lithium  metal=  1  and  a|  i(  v)  =  1  for  pure 


and 


(in  Li(d)  VU(s) 


are  the 


lithium  in  the  solid  state  at  the  electrode  temperature;  and  ff  + ,  ju° 

standard-state  chemical  potentials  of  Li+  ion,  electron,  and  lithium.  The  chemical  potentials  of 
these  species  are  now  expressed  as 


chem 
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Inserting  the  above  information  for  /('*™ne.|ectrolvle)  ,  ,  and  ju^'”  into  Eq.  (131)  leads  to 
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f  RTX 
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(138) 
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Ea'°  =  standard-state  electric  potential  of  the  cell  anode  reaction,  Eq.  (128),  at  the  anode 
temperature,  and  is  given  by 


E 


A,o 


(139) 


Equation  (138)  predicts  that  EAv  would  increase  as  c[  jh  is  increased.  Note  that  the 
mathematical  expression  for  EAv ,  Eq.  (138),  derived  using  the  basic  thermodynamic  principle 
for  the  cell  anode  reaction,  Eq.  (128),  is  similar  to  the  electrokinetics  based  expression  for  EAv 
Eq.  (120),  with  c'J,  =  cu+ .  In  fact,  comparison  of  Eqs.  (138)  and  (120)  leads  to 


(  1  A’1  ^ 
K 

=  exp 

(K>  +K)~Kks) 

kA 

RT 

V  ff  J 

(140) 


Equation  (140)  relates  the  cell  anode  reaction  rate  coefficients  at  a  temperature  and  the 
standard-state  chemical  potentials  of  the  species  involved  in  the  anode  reaction,  //" .  The 
reciprocal  of  Eq.  (140)  leads  to 

KA=e  RT  =  e  RT  (141) 


K 


A,  1 


Notice  that  Eq.  (141)  is  the  mathematical  expression  for  the  thermodynamic  equilibrium  constant 
at  a  temperature,  T[K] ,  where  one  needs  to  know  the  standard-state  Gibbs  free  energy  of  the 
reaction,  Eq.  (128). 


4.1.4  Reversible  or  Open-Circuit  Potential  of  the  Lithium-Based  Cell 

The  cell  reversible  potential  at  a  temperature  is  given  by 

EceU  =EC  -EA 


(142) 


Inserting  the  expressions  for  E^ev  and  EAv  from  Eqs.  (62)  and  (138),  respectively,  into  Eq. 
(142),  leads;  for  the  cell  overall  reaction,  Li  ($)□  Li-in 'the  solid  state  solvent  MyPc'  ,  to 


EceU  = 

rev 
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On  rewriting  Eq.  (143), 
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(144) 

(145) 


(146) 

(147) 


where  AG"  =  M  Pc  i  )  =  standard- state  Gibbs  free  energy  change  for  the  overall  cell 

reaction:  Li  (  v)  [  (Li-intercalated  in  the  solid-state  solvent  MyPc  of  the  cell  cathode) .  Note  that 
as  x  -»  x° ,  yu_in  M  Pc  ->  1 .0  ;  then,  Eq.  (144)  leads  to: 

lim  \_Ecree"~\  =  Ecel,’°  (148) 


This  equation  tells  us  that  with  a  selected  standard- state  of  lithium  in  MyPc,  i.e.,  with  a  selected 

or  defined  value  of  x°,e.g.,  x°  =10“9  g-mole  of  the  intercalated  lithium  per  g-mole  of  initial 

MyPc,  one  can  experimentally  determine  Ecel,’°  at  a  given  cell  temperature,  L] ,  simply  by 

the  measurement  of  voltage  difference  between  the  cell  anode  and  cathode  electrodes  at  the 

cathode  composition  state  of  x° .  Then,  one  should  experimentally  measure  L'f  as  a  function 

of  the  intercalated  lithium,  x ,  in  MyPc  at  the  same  temperature  for  which  EceU'°  has  been 

r  \ 


experimentally  determined.  Using  Eq.  (144)  one  can  compute  yU  m  M  Pc 


as  a  function  of 


f  \ 


One  should,  then,  develop  an  empirical  correlation  for  yLi_m  M  Pc 


or  «u-inMPc  in  terms  of 


\XoJ 


at  one  cell  temperature.  This  procedure  should  be  used  for  at  least  three  temperatures.  The 
empirically  developed  such  correlations  at  different  temperatures  can  then  be  used  to  predict  the 
reversible  cell  voltage  at  any  x  value  during  the  cell  discharge  period.  One  should  be  aware  of 
that  the  cell  reversible  voltage  is  needed  to  predict  the  cell  performance,  i.e.,  actual  power 
provided  by  the  cell  during  its  discharge. 
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4.1.5  Cell  Performance  during  the  Period  of  its  Discharge 


If  the  overall  performance  behavior  of  the  lithium-ion  cell  shown  in  the  sketch.  Fig.  1,  is 
controlled  by  the  lithium  diffusion  in  the  cathode,  MyPc,  material;  then,  the  concentration 
polarization  or  voltage  loss  in  the  electrolyte  can  be  assumed  to  be  negligible  due  to  negligible 
ion  or  electrolyte  concentration  gradient.  At  any  time,  total  cell  current  can  be  predicted  using 
Eq.  (39)  or  (41)  and  total  charge  transferred  from  the  cell  anode  to  cathode  using  Eq.  (40)  or 
(42).  The  lithium  content  in  cathode  initial  material,  MyPc,  during  the  cell  discharge  period  at 
any  time  can  be  computed  from  the  following  result  obtained  using  Eq.  (40)  or  (42) 


v  =  x  + 
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(for  the  cathode  active  material  particle  size  distribution) 


Or, 
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(for  the  cathode  active  material  particles  of  uniform  or  same  size) 

In  Eq.  (149)  or  (150),  xo  =  g-moles  of  the  intercalated  lithium  per  g-mole  of  initial  MyPc  in  the 

cathode  active  material  MyPc  at  the  electrolyte- active  material  particle  interface;  x""  =  initial 
g-moles  of  the  intercalated  lithium  per  g-mole  of  initial  MyPc  in  an  active  material  particle;  and 
N  =  total  number  of  active  material  particle  size  ranges.  One  should  compute  the  cell  reversible 

or  open-circuit  electric  potential  from  Eq.  (144).  The  ohmic  voltage  losses  in  the  cathode 
electrolyte  and  the  cathode  carbon  electronic  conductor  should  be  computed  from  Eqs.  (46)  and 
(48),  respectively.  The  ohmic  voltage  loss  in  the  electrolyte  separator  is  to  be  computed  from  Eq. 
(100).  The  cathode  activation  and  anode  activation  polarizations  or  voltage  losses  should  be 
computed  at  any  time  from  Eq.  (91),  (92),  or  (93)  for  the  cathode  activation  polarization,  ) 

,  and  from  Eq.  (124)  for  the  anode  activation  polarization,  i) sA  .  The  cell  voltage  at  any  time 
during  the  cell  discharge  period  should  then  be  computed  from  the  following  equation 
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(151) 


where  /T'f  =cell  reversible  voltage;  Vre"  —  actual  voltage  difference  between  the  cell  anode 
and  cathode  electrodes;  \^Selectrolyte  |fl ,  | |fl ,  and  |AO“0"  |q  are  the  ohmic  voltage 
losses  in  the  cell  electrolyte  separator,  cathode  electrode  electrolyte,  and  cathode  carbon 
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electronic  conductor,  respectively;  and  |(-//f  )|  and  |//;'  |  are  the  cathode  and  anode  electrode 

activation  polarization  voltage  losses,  respectively.  The  actual  electric  power  delivered  by  the 
cell  during  the  cell  discharge  period  at  any  time  is  given  by 

p  _  jt-ycell  _  jt,C-ycell  (152) 

for  the  cell  current  I‘  =  I'  c  (current  limited  by  the  lithium  diffusion  in  the  cathode  active 
material  particles).  The  thermodynamic  efficiency  at  any  time  may  be  defined  as 


jt  ,C-y  cell  -ycell 

7  thermody  ~  p,C  pcell  ~  pcell 


(153) 


4.1.6  Cathode  Active  Material  Spherical  Particle  Utilization  to  Store  Lithium  (i.e. 
Coulomb  Efficiency  of  a  Cathode  Active  Material  Particle)  during  the  Cell 
Discharge  Period 

The  formulation  to  predict  the  transient  lithium  concentration  profiles  in  a  spherical  particle  of  a 
cathode  active  material  as  a  function  of  the  dimensionless  radial  distance  is  provided  for  a 
constant  lithium  ion  flux  or  a  constant  electrolyte-cathode  active  material  particle  interfacial 
current  density.  This  formulation  is  then  employed  to  develop  the  mathematical  expression  to 
predict  the  coulomb  efficiency  of  the  cell  cathode  consisted  of  uniform- sized  active  material 
spherical  particles  under  the  condition  of  lithium  diffusion  in  the  particles  controlling  the  overall 
cell  discharge  process. 

The  partial  differential  equation  (PDE)  describing  lithium  diffusion  in  a  spherical  particle  of  the 
cathode  active  material  is 


dc  D  d  f  dc^ 
dt  r 2  dr  \  dr , 

The  initial  and  boundary  conditions  are: 

t<0,c  =  cini 
t  >  0,  r  -  0,  c  =  finite 

or,  the  boundary  condition  represented  by  Eq.  (156)  can  also  be  given  by 


(154) 


(155) 

(156) 


t  >0, 


^  dc 

\ 

r=oJ 

=  0 


(157) 


(symmetry  of  the  intercalated  lithium  concentration  with  respect  to  the  radial  distance  at  the 
spherical  particle  center) 
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The  second  boundary  condition,  at  r  =  a  (particle  surface),  is  given  by 


4C 

=  Js=^r  (158) 

.  r=a)  F 

where  D  =  lithium  mass  diffusivity  averaged  over  lithium  concentration  in  the  cathode  active 
material  particle,  Js  -  lithium  ion  molar  flux  at  the  particle  surface;  assumed  here  to  be  constant, 
|/sc|  =  magnitude  of  the  electrolyte- active  material  particle  interfacial  current  density;  assumed  to 
be  constant,  F  -  Faraday’s  constant  =  96487  coulomb  per  g-mole  of  Li+  (or  per  g-equivalent). 

The  solution  to  the  PDE,  Eq.  (154),  for  the  initial  and  boundary  conditions,  Eqs.  (155)  through 
(158)  is  adopted  from  [3]  and  is  given  below. 
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Ut  v 

where  r  =  —  ,  g  =  —  .  x  =  moles  of  the  intercalated  lithium  per  mole  of  the  cathode  active 
a ~  a 

material  at  any  location  E,  in  the  active  material  particle,  xini  =  initial  moles  of  lithium  per  mole 
of  the  cathode  active  material,  pcathode  =  cathode  active  material  particle  density,  and  Mcmhode  = 
mol.  wt.  of  the  cathode  active  material.  The  ( cian )  s  are  the  roots  of  the  following  mathematical 
expression: 

(aan)cot(aan)  =  l  (160) 


Positive  roots  of  Eq.  (160)  are 

\aa{)  =  4.4934,  (aa2)  =  7.7253,  (aa3)  =  10.9041, 
(<7 a4  )  =  14.0662,  (aa5 )  =  17.2208,  etc. 


(161) 


At  the  particle  center,  i.e.,  at  %  -  0 ,  Eq.  (159)  reduces  to: 


D(C center  Cini ) 


J  .a 


lHPenu,-Xmi)p< 


cathode 


(J,a)M 


cathode 


n= oo  exp  l—T(aan)2) 

3r_o  3-2V _ 1 _ - 

n= i  (fl«„)sin(aajI) 


(162) 


Equation  (162)  predicts  the  lithium  concentration  at  the  spherical  particle  center  as  a  function 
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f  Dt 

or  r  =  — r  . 


At  the  particle  surface,  i.e.,  at  g  =  1 ,  Eq.  (159)  leads  to: 


Die  —c) 

\  s  ini  ) 


J  a 


_  ^  ^ ini  )  P cathode 


cathode 


3r  +  0.2-2^ 


exp  j 


(163) 


Equation  (163)  predicts  the  lithium  concentration  at  the  spherical  particle  surface  as  a  function  of 

f  Dt) 

the  dimensionless  time,  r  =  —  . 

I  a  J 

If  the  time  for  lithium  concentration  at  the  cathode  active  material  particle  surface  to  reach  cssa,d 
is  tssa,d  ,  then  the  amount  of  lithium  that  would  have  accumulated  in  the  particle  by  this  time  is 
given  by 


n: 


t -actual 
Li-surface  sat'd 


=  Ana  Jt 


sat'd 


(164) 


It  should  be  noted  that  when  cs  becomes  equal  to  cssatd  during  the  cell  discharge,  cell  voltage 

would  drop  steeply  and  the  cell  then  must  be  declared  discharged.  The  amount  of  lithium 
required  to  be  transferred  or  intercalated  to  uniformly  saturate  the  spherical  particle  is  given  by 


N, 


t  -sat'd 
Li 


f4  3^ 

(xsat'd  -JC  ' 

\  P cathode-active  material 

v3  ) 

M 

cathode-active  material  _ 

(165) 


The  cathode  active  material  utilization  or  its  coulomb  efficiency  to  store  lithium  during  the  cell 
discharge  period  is  defined  as 


U  = 


N, 


t  -actual 

Li-surface  sat'd 


N, 


t  -sat'd 
Li 


Using  the  information  from  Eqs.  (164)  and  (165)  in  Eq.  (166)  leads  to: 


U  =3 


J. 


f  j.  sat'd  \ 


v  a  J 


(sat'd  \ 

^  ^ ini  J  /^cathode- 

M 


•active  material 
cathode-active  material 


(166) 


(167) 
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Equation  (167)  is  now  expressed  as 


U  = 


3r 


sat'd 


y  (Jsa)M  cathode  J 


(168) 


From  Eq.  (163), 


3  zfd  +0.2- 

n= OO  exp| 

sat'd  (  _  \ 2 

(aa„) 

I 

v  (Tia)Mcath<xle  ^ 

n= 1 

Inserting  the  information  from  Eq.  (169)  into  Eq.  (168)  leads  to 


o  _  sat'd 

3r.s 

exp  (-rf"y  (aan  f  j 

■< 

3zsa,d  +0.2  2Y  1  ,  ' 

> 

»= 1  \aan) 

(169) 


(170) 


where,  r  = 


Dt 


sat'd 


Numerical  assessment  of  ^exp 


_  sat'd 


(aan) 


2  A 


(aa„y 


indicates  that  this  sum  makes  negligible 


contribution  to  the  function  shown  in  the  denominator  of  the  right-hand  side  of  Eq.  (170)  for 


Tssatd  >  1 .  So,  Eq.  (170)  is  written  as 


U  = 


3  T 


sat'd 


3  Tsa,'d  +0.2 


(171) 


valid  for.  zs“,d  >  1 . 

One  should  use  Eq.  (170)  to  compute  U  for  rssatd  <1  and  Eq.  (171)  for  zssatd  >  1 . 

Once  U,  the  cathode  active  material  particle  collection  efficiency,  has  been  determined  using  Eq. 
(170)  or  (171)  for  a  fixed  tssa,d  and  a  particle  radius;  then  one  should  compute  the  lithium  ion 
flux  at  the  electrolyte-active  material  interface,  Js ,  from  Eq.  (168).  Then,  corresponding  to  the 
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computed  J  s  value,  one  should  determine  the  electrolyte-cathode  active  material  particle 
interfacial  current  density,  /sc . 


C  =  J.F 


(172) 


To  show  the  effect  of  the  cathode  active  material  particle  radius  on  U,  Js ,  and  ,  one  should 
compute  them  for  a  selected  fa,d  value  (e.g.  tssa,d  =  3600  seconds )  for  a  given  cathode  material 
as  a  function  of  particle  radius,  a  or  r  d  (dimensionless  time  for  the  particle  surface 
saturation).  For  the  practical  application  of  such  computed  data  in  the  design  and  performance 
analysis  of  a  cathode  of  a  given  material,  one  should  develop  plots  of  U,  Js ,  and  ics  as  a  function 

of  Tssa,d  and  the  active  material  particle  diameter,  2a,  for  ts‘"d  =  3600  secs . 


In  what  follows  the  development  of  the  formulation  to  compute  the  transient  lithium 
concentration  profiles  in  a  cathode  active  material,  spherical  particle  for  the  short  time  condition 
is  briefly  presented.  This  formulation  should  only  be  used  when  the  cell  discharge  period  is  so 
short  that  the  intercalated  lithium  concentration  at  the  particle  surface  would  have  reached  the 
saturation  level  while  the  lithium  penetration  toward  the  particle  center,  via  the  molecular 
diffusion  process,  is  still  in  the  spatial  region  very  near  the  particle  surface. 

The  flux  distribution  of  the  intercalated  lithium  inside  the  cathode  active  material  particle  is 
given  by  the  following  partial  differential  equation  (PDE)  [8]. 


*=-±i.(rv) 

dt  r°  dr '  r 


(173) 


where  c  =  molar  concentration  of  the  intercalated  lithium  at  any  radial  position,  r,  at  any  time,  t, 
inside  the  cathode  spherical  particle  and  J  r  -  lithium  diffusion  flux  along  the  radial  position;  it 
is  positive  if  lithium  diffuses  away  from  the  particle  center  and  negative  if  its  diffusion  is  toward 
the  particle  center.  The  radial  diffusion  flux  J  r  is  expressed  as 


J  =  D 


dc 

dr 


(174) 


Differentiating  the  PDE,  Eq.  (173),  with  respect  to  r  and  using  the  expression  for  Jr  in  Eq.  (174) 
with  constant  D  leads  to 


=  D  — 

dt  dr 


1  5  (  2 
—  —  r 
r2  dr v 


Jr) 


(175) 


For  the  situation  of  short  times,  lithium  ions  after  reduction  and  entering  the  active  material 
particle  diffuse  (as  Li+  +e~)  toward  the  particle  center.  It  is,  however,  assumed  that  they  stay 
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near  the  particle  surface  until  the  lithium  concentration  at  the  particle  surface  located  at  r  =  a 
reaches  the  saturation  level  of  cssa,d  or  xssatd ;  where  cssa,d  =  saturated  molar  concentration  and 
xs'u d  =  g-moles  of  lithium  per  g-mole  of  the  cathode  active  material  at  the  saturation  level.  The 
intercalated  lithium  ((Li+  +  e  j  =  Li  j  atoms  ‘see’  the  particle  spatial  region  away  from  “a  very 

thin  spatial  region  near  the  particle  surface”  as  if  it  is  not  approachable  during  the  cell  discharge 
period  until  the  lithium  concentration  at  the  particle  periphery  reaches  the  saturation  level.  For 
this  situation,  the  lithium  diffusion  inside  the  cathode  active  material  particle  is  very  near  the 
particle  surface.  Therefore,  the  PDE,  Eq.  (175)  is  rewritten,  with  the  assumption  of  r  — >  a,  as 
follows. 


dt  dr 


1  d  (  2 1  \ 

—  ^{a  J> 
a  or  v  ' 


(176) 


Equation  (176)  becomes 


dt  dr 


—  Jr 

dr 


(177) 


Position  of  a  spatial  point  inside  the  cathode  active  material  particle  is  expressed  as 


y  =  a-  r 


(178) 


where  y  —  distance  of  the  spatial  point  measured  from  the  particle  surface  and  r  =  distance  of 
the  same  spatial  point  measured  from  the  particle  center  along  a  particle  radius.  Note  that  at  r  = 

a,y  =  0  and  at  r  =  0,  y  =  a  and  —  =  -1 . 

dr 


Now, 


dr 


'  dJr  dy  \  _  dJr 
y8y  dr  J  dy 


knowing, 


5 

_  d 

d 

’  5 

' 8Jr ' 

dy 

dr 

dr  _ 

dr 

dy  _ 

dr 

_^y  _ 

_^y 

k  dy  j 

dr 

Hence, 


dr 


’  5 

"a/,," 

ay 

dy 2 


(179) 


(180) 
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Inserting  the  information  from  Eq.  (180)  for 


into  Eq.  (177)  leads  to 


dr2 


dJ_ , 

dt 


=  D 


d_% 

dy2 


(181) 


Note  that  J r  (  =  \j r  |)  has  a  positive  number  value  with  respect  to  y. 

|4C| 

Divide  both  sides  of  Eq.  (181)  by  Js  =  J — 1 ,  the  lithium  molar  flux  at  the  particle  surface  along 

F 

the  radial  distance  toward  the  particle  center  corresponding  to  a  fixed  electrolyte- active  material 
particle  interfacial  current  density  I ics  I .  The  PDE,  Eq.  (181)  becomes 


d2r 

w 


(182) 


where 


r  = 


(183) 


The  initial  and  boundary  conditions  to  solve  the  PDE,  Eq.  (182),  are 
I.C.:  For  t<  0, 


r  =  0,  for,  0  <  y  <  a 

(184) 

For  t>  0, 

r  =  l,  aty  =  0  (at  the  particle  surface) 

(185) 

For  t>  0, 

r  =  0,  aty  — >"qo" 

(186) 

Here  the  PDE,  Eq.  (182),  is  solved  using  the  method  of  combination  of  independent  variables.  A 
dimensionless  quantity,  i) ,  is  defined  as 

(187) 

Now  the  PDE,  Eq.  (182),  and  the  initial  and  boundary  conditions  are  transformed  into  the 
dimensionless  form. 
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ar_jr  d  y  j  _  dr  [  y  \d  f  1 

dt  dr]  dt  dr/  6t  {  y4/;/  )  dr]  \  ^4/)  jdt{sj't 

-  '*r)(  y  ¥  O-f  7Vr 

2 1)  \  2t)dr] 


(188) 


<3T  _  dT  dr]_  _dT_  d_(  y 
dy  dr]  dy  dr]  dy{^jjf)i 


dr  1 


drj{y/4Dt  J\dy)  {j4Dt)dr] 


1  )dr 


(189) 


Differentiating  both  sides  of  Eq.  (189)  with  respect  to  y, 


e2r  d 


1  )dr 


1  \d(dT 


dy2  dy  \[ j4Dt J  dr]  \  y  ^ADt  j  dy  {dr/ 
1  \d(dT^dr]~ 

mUvUv J  dy  _ 


thus, 


dT  _  1  d T  1  _(  1  WT 
dy2  y/4, Wt  dr f  ^/4Dt  U Dt  J  dr] 


(190) 


(191) 


dr  ^2  j— < 1 

Inserting  the  mathematical  expressions  for  —  and  — j  from  Eqs.  (189)  and  (190), 

dt  dy 

respectively,  into  Eq.  (182)  and  simplifying  leads  to 


d  r  „  dr 


-t+2,i—  =  o 


dr/  dr] 


(192) 


The  initial  and  boundary  conditions  in  the  transformed  form  are 


B.C.l: 


B.C.2: 


at  77  =  00,  r  =  0 


at  77  =  0,  r  =  1 


at  77  =  00,  r  =  0 


(193) 


(194) 


(195) 


Noting  that  the  initial  and  B.C.2  conditions  are  identical,  the  above  conditions  are  now  rewritten 
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B.C.l: 

B.C.2: 


at  tj  =  0,  r  =  1 
at  7j  =  go,  r  =  0 


(196) 

(197) 


Now,  the  differential  Eq.  (192)  is  solved  using  the  boundary  conditions,  Eqs.  (196)  and  (197). 
Let, 


dT 

dr/ 


=  ¥ 


(198) 


Therefore, 


d2 T  _  dy/ 
dr/2  dr/ 


(199) 


Inserting  the  information  from  Eqs.  (198)  and  (199)  for  — 

dt] 

(192)  leads  to 

and  — ^  ,  respectively,  into  Eq. 
dr/ 

^  +  2W  =  0 
dr/ 

(200) 

Separating  the  variables, 

— -  =  -2  r/dr/ 

¥ 

(201) 

Integration  of  both  sides  of  Eq.  (201)  leads  to 

l°gf  ¥  =  ~rJ2  +  A 

(202) 

Simplification  of  Eq.  (202)  leads  to 


y/  =  e~’l2+A  =(eA')e-n2  =  Ae^2 


(203) 


Inserting  the  expression  for  (//from  Eq.  (203)  into  Eq.  (198)  leads  to 


dT 

dr/ 


Ae 
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(204) 


(205) 


dF  =  Ae  n  dr/ 


Integrating  both  sides  of  Eq.  (205), 


V 

T  =  AjV ^dri  +  B 


(206) 


Rewriting  Eq.  (206), 


w-// 

Y  =  A  |  el‘du  +  B 


(207) 


where  u  =  ’dummy  variable’  for  r/ ,  and  B  =  second  constant  of  integration.  Now,  the  integration 
constants  A  and  B  are  evaluated  using  the  boundary  conditions,  Eqs.  (196)  and  (197). 


U=7] 

(rl  )  =  A  f  e 

\  I  at  u=rj= 0  )  J 


dull  +B 


(208) 


I  at  u=ij= oo 


u=i 

H  1 


e  "  du  +B 


(209) 


The  expressions  in  Eqs.  (208)  and  (209)  lead  to 


1  =  A  |  e~u  du  +  B 


(210) 


U  —  JJ 

0  =  A  j  e  u  du  +  B 


(211) 


The  solution  to  Eqs.  (210)  and  (211)  for  A  and  B  leads  to 


(212) 


A  =  — 


|  e  "  du 


(213) 
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Substituting  for  B  and  A  from  Eqs.  (212,  213)  into  Eq.  (207)  leads  to 


r  =  i— 


|  e~u  du 
_o _ 

oo 

|  e~u  du 


(214) 


2 


=  1- 


4n 


|  e  du 


f  2 

ysfx  o 


f—  f  e  "  du 


(215) 


because  erf  ( oo)  =  l. 


Rewriting  Eq.  (216), 


r=i- 


f  e  du 

I  o 


sfx  J  erf  (oo) 


2  n 

=  1 - ^=je du 


r  =  y1  =  1  ~erf(Tj)  =  erfc(rj) 


(216) 


(217) 


2  n 

where  the  error  function  of  77,  erf  (77)  =  -=  I  e  du  and  the  complementary  error  function  of 

dX  i 

1 1 . 

2  M=CO 

erfc(r/)  =  —~=  |  e~u  du  =  \- erf  (ij )  and 

u=v 


77: 


y 


^/4D7 


(y/g)  _  (.v/fl)  _  y 

74  15 

V«7  V  2 


!D7 
fl2 


(218) 


y 

where  the  dimensionless  distance  Y  =  —  and  r  = 

a 


(  Dt  ' 


v«  y 


,  dimensionless  time. 
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One  can  predict  T 


f  J > 

r  Dt ' 

_  r 

as  a  function  of  /;  or  Y  at  each  r  = 

V  Jsj 

l  a~  J 

or  t  value  using  Eq.  (217). 


At  77  =  2 ,  erfc(rj)  =  0.01 .  Therefore,  a  mass  transfer  boundary-layer  thickness,  8m ,  is  defined 
corresponding  to: 


n 


Y  _  y  5rn  _2 

2\fr  2<i \fr  2a  4z 


The  mass-transfer  boundary-layer  thickness,  Sm ,  is  then  given  by 

8m  =  Aa4r  =  4 


(219) 


(220) 


The  mass-transfer  boundary-layer  thickness,  8m ,  is  a  natural  length  scale  for  the  lithium  ion 
mass  flux  in  a  cathode  active  material,  spherical  particle.  The  distance,  8m ,  is  a  measure  of  the 
extent  to  which  lithium  ions  have  penetrated  into  a  particle  of  the  cell  cathode  active  material  in 
the  elapsed  time,  t.  It  is  proportional  to  \ft  . 


To  obtain  the  lithium  concentration,  c,  distribution  in  a  spherical  particle  of  the  cathode  active 
material,  the  lithium  ion  flux  inside  the  particle  toward  the  particle  center  is  expressed  by 


Jr=-D— 

dr 


Because,  r-a-y;  therefore, 


dr  =  —dy 


(221) 


(222) 


From  Eq.  (221), 


dc 

(-dy) 


Rearranging  Eq.  (223)  and  integrating  the  resulting  expression  leads  to 


(223) 


(224) 


(225) 
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Or, 


(226) 

(227) 

(228) 

(229) 

(230) 

(231) 

(232) 

(233) 
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(234) 


=  J, 


16/  V 


1 

—  oo 


nD 


I  I 


e  "  du 


di] 


=  J. 


16/ 

nD  , 


(  =0° 

1 

C=T) 


j"  e  du 


dC 


=  J 


^  16/  V 

u=  00  I 

f 

\(r  >  1 

f  e  z/C 

[nDj 

J 

U=J]  \ 

[A  J 

=  J. 


16/ 
nD . 


M=00 

I 


du 


=  J, 


16/ 

nD 


u=  oo 

J  jV"  (z/  —  77)| 


u=rj 


=  J. 


16/ 

nD . 


w=co 

J  {Me_ 


-77c  "  >z/zz 


Therefore, 


{C~Cini)  =  Js 


16/ 

nD 


il  —  JJ 

J  we"  du  - 1]  J 


du 


with  u  =  z,  zzz/zz  = 


z/z 


00 

r  m2  7 

[VO 

1  ue  du  = 

7 

l  2  J 

(235) 

(236) 

(237) 

(238) 

(239) 

(240) 

(241) 

(242) 


Substitution  for 


from  Eq.  (242)  into  Eq.  (240)  leads  to 
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(243) 


(244) 

(245) 

(246) 

(247) 

(248) 

(249) 
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where  (the  complementary  error  function  of 


f  Y  ) 

f  Y  > 

f  2  3  “7 

=  erfc 

(2  s/t  j 

itf]  l 

2-Jt 


The  equation  (249)  is  now  written  as 


(' C~Cini)D 

(  2  ^ 

sfr 

ryri 

Y  erfc 

f  Y  } 

Jsa 

vv/^y 

f  y2) 

pvn  - 

2  J 

[2yfr  y 

v4ry 


(250) 


(251) 


(valid  for  r>  0;  0<  Y  <0.1  to  satisfy  the  condition  of  short  distance  penetration  by  lithium,  i.e. 
Li++e"  pair,  from  the  electrolyte-active  material  interfacial  surface  into  the  particle). 


At  the  particle  surface,  i.e.,  y  =  0  or  F  =  —  =  0 ,  c  =  cs,  Eq.  (251)  becomes 

a 


(c,  -cm)D 

r  2  y 

Jsa 

vv/^r  J  _ 

\fr 

T 


2^(0)er/c(0) 


(252) 


(—) 

y 


y 


sfr 


(253) 


1 


Jr 


(254) 


Notice  that  the  particle  surface  concentration,  i.e.,  (c5  -  cini )  is  proportional  to  \fr  .  Expressing 
the  equations  (251)  and  (254)  as  follows, 


(C~Cuu)D 

x ini  ^  p cathoc[eD 

f  2  ] 

Jsa 

M  cathode  (j  sa) 

1 4n) 

v/r 


fVTl 


Y  erfc 


f  Y  ) 

v2v/r  y 


(255) 


(c  —c  )D 

\  s  ini  J 

( ^ ini  )  P cathode  ^ 

(  2  ^ 

Jsa 

M cathode  {Jfi) 

1 J 

v/F 


(256) 
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At  the  time,  t  =  tssa,d  ,  when  the  lithium  concentration  at  the  particle  surface,  r  =  a  or  Y- 1 
reaches  the  saturation  level  for  a  given  cathode  active  material  of  density,  pcathode ,  and  molecular 


weight,  Mcathode ;  then,  c  =  cfd  and  xs  =  xfd  and  r  =  r  fd  = 


Dt 


sat'd 


a 


Equation  (256)  is  rewritten  for  the  cathode  active  material  surface  saturation  with  lithium. 


(csa,'d -c.  )D 

\  s  ini  J 


J  a 


(Xf<l  ~Xini)PcathodeD 


M . 


( J  A  ) 


f  o  A 


j 


I  sat'd 


(257) 


Note  that  in  Eqs.  (255)  and  (256),  x  =  g-moles  of  Li  per  g-mole  of  cathode  active  material. 


Equation  (255)  can  be  used  to  predict 


at  each  r  and  Eq.  (256)  should  be  used  to  predict 


(C~Cini)D 

ini  )  Pcathode ^ 

Jsa 

or 

M  cathode  (J  sa) 

(Cs~Cini)D 

(■*)  Xini  )  PcathodeD 

J  a 

M cathode 

as  a  function  of  Y 


as  a 


function  of  t.  For  a  fixed  value  of  tssa,d ,  e.g.,  tssa,d  -  0.001  second  for  the  lithium  concentration 
at  the  particle  surface  to  reach  [xssMd  -  xM  j  ,  one  can  compute  ( Jsa )  from  Eq.  (257).  Then,  one 


can  compute  Js  corresponding  to  a  value  of  the  spherical  active  material  particle  radius,  a.  One 
may  call  this  lithium-ion  flux  at  the  electrolyte -active  material  interface  Jssa,d  ,  i.e.,  J ,  =  Jssa,d  , 
lithium  ion  flux  required  at  the  particle  surface  for  the  lithium  concentration  at  the  particle 
surface  to  reach  the  saturated  level  in  a  desired  period  of  time,  tssa,d  or  rs'""f .  One  can  then 
compute  the  electrolyte-active  material  interfacial  current  density  to  satisfy  this  requirement,  i.e., 
|/f  |  =  j°at'dF  for  a  selected  cathode  active  material  particle  radius.  It  is,  here,  reminded  that  the 

cell  should  be  considered  discharged  when  the  lithium  concentration  inside  the  particle  at  its 
surface  reaches  the  saturation  level. 


The  cathode  coulomb  efficiency,  U,  to  store  lithium  in  a  particle  of  a  cathode  active  material, 
(see  Eq.  168),  is  given  below. 


U  = 


3r 


sat'd 


D(xT'd -x„,)p,M 


(J,a)M 


cathode 


(258) 


Inserting  the  expression  for 


D(xfd  -xini)pcalhode} 

(Jsa)Mcathode 


from  Eq.  (257)  into  Eq.  (258)  leads  to 
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(259) 


U  = 


V  2  J 


sat'd 


sat'd 


So,  one  can  compute  the  cathode  coulomb  efficiency,  U,  from  Eq.  (259)  for  a  zs‘"d  =  —h, —  for 


a 


a  fixed  tftd  value  for  the  various  values  of  the  cathode  active  material  particle  radius,  a,  to 
observe  the  effect  of  the  particle  dimension  on  the  coulomb  efficiency  by  the  use  of  Eq.  (259) 
derived  for  the  short-time  cell  discharge  period. 


Properties  of  erf{rj )  and  erfc 
for  their  use  in  Eq.  (255). 


n 


2yfr 


;  and  some  of  their  tabulated  values  are  given  below 


r\  1  f  rs  \  °° 

erf  (7)  =  -7=  f  e  “  du,  erfc  (rj)  =  l- erf  (rj)=  -=\\e~^du 

dxJ0  W  xJi 


erfM  =  ^n-  — 
erfc(rj)  =  1  - 


3  5 

77  +  77 


V 


5-2!  7-3! 


4ft 


V 


if  77 5 

77  -  +  — - 

3-1!  5-2!  7-3! 


V 


erf  (-77)  =  ~eif  (77) ,  <?t/(0)  =  0,  ct/  (00)  =  1;  erfc  {ft)  =  1,  er/c(oo)  =  0. 


Also,  erf  (77)  □  1  — 


erfc  [if)  □ 


shin 


x _ 1  |  1-3  1-3-5  | 

V  (2?72)2  (2t72)3 


f 


'[nil 


x _ 1  (  1-3  1-3-5  | 

V  (2?72)2  (2t72)3 


V 


Table  1.  Some  Values  for  the  Error  and  Complementary  Error  Functions 


Y 

71  ~  2  4r 

erf  (77) 

erfc  {if) 

0.00 

0.00000 

1.000000 

0.05 

0.056372 

0.943628 

0.10 

0.112463 

0.887537 

0.15 

0.167996 

0.832004 
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0.20 

0.222703 

0.777297 

0.25 

0.276326 

0.723674 

0.30 

0.328627 

0.671373 

0.35 

0.379382 

0.620618 

0.40 

0.428392 

0.571608 

0.45 

0.475482 

0.524518 

0.50 

0.5205 

0.479500 

0.55 

0.563323 

0.436677 

0.60 

0.603856 

0.396144 

0.65 

0.642029 

0.357971 

0.70 

0.677801 

0.322199 

0.75 

0.711156 

0.288844 

0.80 

0.742101 

0.257899 

0.85 

0.770668 

0.229332 

0.90 

0.796908 

0.203092 

0.95 

0.820891 

0.179109 

1.00 

0.842701 

0.157299 

1.10 

0.880205 

0.119795 

1.20 

0.910314 

0.089686 

1.30 

0.934008 

0.065992 

1.40 

0.952285 

0.047715 

1.50 

0.966105 

0.033895 

1.60 

0.976348 

0.023652 

1.70 

0.98379 

0.016210 

1.80 

0.989091 

0.010909 

1.90 

0.99279 

0.00721 

2.00 

0.995322 

0.004678 

2.10 

0.997021 

0.002979 

2.20 

0.998137 

0.001863 

2.30 

0.998857 

0.001143 

2.40 

0.999311 

0.000689 

2.50 

0.999593 

0.000407 

2.60 

0.999764 

0.000236 

2.70 

0.999866 

0.000134 

2.80 

0.999925 

7.5E-05 

2.90 

0.999959 

4.1E-05 

3.00 

0.999978 

2.2E-05 

4.1.7  Determination  of  Voltage  Loss  Associated  with  the  Diffusion  Transport  of  Lithium, 
i.e.,  the  Lithium  Diffusion  Voltage  Polarization  in  the  Cathode  Active  Material 

Diffusivity  of  lithium  ion  in  the  cathode  solid-state  active  material  as  the  solvent  for  the  solute, 
lithium  ion  is  related  to  its  ionic  conductivity,  ,  via  the  following  expression  [10]: 
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(260) 


D  +  = 

Li 


RTkc 


h-F)r 

(the  Nemst-Einstein  relation) 


Subtracting  Eq.  (162)  from  Eq.  (163)  leads  to: 


(Cs  C  center}  D 


J  a 


Os  ~Xce„ter)D 


J  a 


P cathode 
V  M  cathode  J 


0.5 +  2X 


exp(-r(aan  j  j  exp 


(aan)sm(aan)  ( aanf 


(261) 


Substitution  for  D  (  =  D  )  from  Eq.  (260)  into  Eq.  (261)  and  simplification,  with  Js  =  — , 

Li  V  /  F 

leads  to 


c-c 


center 

c  j 


RT 


V  r  J 


•  C 

5 

h 

a 

V 

4 

-  ) 

x—x 


center 

x  j 


^ (  RT' 

FF; 


•  C 

\ 

h 

Cl 

< 

V 

+  J 

0.5  +  2X 


n- 1 


exp  -r(aan^J  J  exp  -r(acfn) 


(aan  )sin  (aan )  (aa„)2 


(262) 


where  the  lithium  concentration  j  averaged  over  the  dimensionless  particle  radius  at  any  time 
x  is  given  by 


1.0 

x  =  |  xcU;  (263) 

0 

Note  that  v,  at  any  time  t,  is  to  be  determined  as  a  function  of  \  using  Eq.  (159).  Then,  Eq.  (263) 
is  to  be  employed  to  obtain  x  at  each  x. 
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Defining 


T- 

—  X 

center 

' RT v 

y c 

C  center 

X 

r rtX 

_v 

X 

J 

_V 

c 

J 

=  (AVd)C,  voltage  loss  associated  with 


the  lithium  ion  diffusion  inside  a  cathode  active  material  spherical  particle  due  to  the  situation  of 
xs  and  xcenter  being  different  from  x  ;  then,  Eq.  (262)  is  written  as 


f 

•c 

X 

l 

Cl 

s 

k. 

v  +  J 


0.5  +  2  X 


n= 1 


exp(-r(<m„)  j  exp(-r(<mH) 


(aan)sin(aan)  (aa„)2 


(264) 


where  |/sc|  =  constant  electrolyte-cathode  active  material  interfacial  surface  current  density  during 

kc  = 

the  cell  discharge  period,  a  “cathode  material  particle  radius,  and  +  lithium  ion 

(a  yD)c 


conductivity  in  the  cathode  active  material.  Equation  (264)  can  be  used  to  determine 
voltage  loss  associated  with  the  lithium  ion  diffusion  in  the  cathode  active  material  as  a  function 
of  r  during  the  cell  discharge  period.  For  the  general  application  of  the  result  in  Eq.  (264)  for  the 
spherical  particles  of  a  cathode  active  material,  it  is  suggested  that  the  dimensionless 

lithium-diffusion-voltage  loss,  |  f  \if  |  a  ]  | ,  should  be  calculated  as  a  function  of  r  and  a  plot  of 

(AVj 


K+ 

v  +  J 


f 

. c 

X 

l 

a 

s 

V  +  J 


vs.  t  be  developed. 


4.1.8  Evaluation  of  the  Effect  of  Lithium  Storage  on  the  Cell  Reversible  Voltage 

The  cell  reversible  voltage  is  given  by  Eq.  (144).  If  the  standard- state  lithium  composition,  x° , 
is  unknown,  one  may  eliminate  it  from  the  cell  voltage  expression  by  first  writing  Eq.  (144)  for  a 
known  value  of  E^\  at  a  lithium  composition,  x  =  x1  value  and  then  subtracting  the 

expression  for  from  Eq.  (144)  to  obtain,  on  simplification, 


pi  cell  pi  cell 

rev  ^rev,  1 

RT^ 


=  -In 


Yu 

Yua  ) 


rx A 


vxiy 


(265) 


Note:  The  value  of  the  function  within  [  ]  should  be  >  1. 
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Notice  that  the  reversible  cell  voltage  loss  can  be  determined  as  a  function  of  x  provided 


f  \ 
Yu 


is  known  as  a  function  of  x  or 


\X\J 


Uu  ij 

during  the  cell  discharge  period.  To  evaluate  the  effect  of 


lithium  storage  on  the  reversible  voltage,  x  =  xs  (at  the  particle  surface);  therefore,  Eq.  (265) 
becomes 


■pi  cell  pi  cell 

^  rev  ^rev,  1 


( RT\ 


V  r  J 


=  -In 


r  \ 

Yu 


Yu,  i ) 


fx' 


V*i7 


=  F(x.) 


(266) 


xs  should  be  computed  from  the  following  equation  derived  from  Eq.  (163). 


T  =  Xini  + 


(J,a)M, 


cathode 


Dp 


cathode 


n= oo  exp  l—T(aan)2) 

3r  +  0.2-2X— ( - 75  ~ 

»=1  P1®,,) 


(267) 


Recognizing  Js  =  — ,  Eq.  (267)  is  now  expressed  as 
F 


Xs  =  Xi„i  + 


F 

v  J 


M 


cathode 


Dp 


cathode 


3r  +  0.2-2^ 


exp 


(268) 


If  the  correlation  for  yu  ,  activity  coefficient  of  Li  (in  fact,  Li+),  as  a  function  of  the  intercalated 
or  moveable  lithium  ion  composition  is  not  available;  then,  F (xs )  at  an  xs  =  x1  +  Sx  can  be 
obtained  from  the  following  derived  Taylor  series,  provided  Sx  has  such  value  that  yu  (at  an 
xs  =  x,  +  Sx )  does  not  differ  from  yu  j  at  x,  to  any  significance. 


FU) 


' Sx ' 

rn 

'  Sx^ 

2 

'  8x^ 

3 

'  Sx' 

vh  J 

vXi  V 

+  | 

v3!y 

vxi  y 

^  4!y 

v*i  V 

(269) 


where  the  remainder  R5 


is  given  by 


"24^ 

v5!v 

U  J 

(the  Lagrange’s  form) 


(270) 
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5x 


(271) 


R*  = 


{(xl+sx)-^y 


(the  Cauchy’s  form) 


with  Xj  <  g  <  (xj  +  Sx)  . 


To  observe  the  effect  of  lithium  storage  in  the  cathode  active  material  particles  on  the  reversible 
voltage,  F(xs)  should  be  computed  as  a  function  of  x  using  Eq.  (268);  then,  the  plot  of 


ricell  pi  cell 
^  rev  rev, l 

~rRT^~ 


vs.  x  should  be  developed. 


Equation  (256)  developed  to  predict  the  lithium  surface  concentration  in  a  cathode  active 
material  for  the  condition  of  lithium  penetration  in  a  short  distance  from  the  electrolyte-active 
material  interface  may  also  be  used,  in  approximation,  to  compute  the  voltage  loss  associated 
with  the  lithium  diffusion  inside  the  cathode  active  material  particles  and  to  observe  the  effect  of 
lithium  storage  on  the  cell  reversible  voltage. 

The  lithium  ion  diffusivity  D  in  Eq.  (256)  is  expressed  in  terms  of  the  lithium  ion  conductivity 
using  the  Nernst-Einstein  relation,  Eq.  (260)  to  obtain: 


(c* 

-c„)RT«rf 

|'jj! 

_  ini  )  P cathode 

RTk\ 

1 

(vf) 

1 F 

Jsaj 

M  cathode  (Jfl)  1 

(7>) 

1 F 

V?  =  1.1284 >/r 


Because  cT^+  =  x, :+  Pcathode  ;  so, 


Li  Li  ,avg  ^ 
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(  ^ ini  )  P cathodi 

M  cadu>de(J >) 


RTkc 


Li  ,avg  m 


Pcathode  {F  •  Fj 


=  1.1284 yfr 


cathode 


Xs~Xini 


V  Xavs  J 


f  RT  \  kc 


V  r  J 


F{Jsa ) 


=  1.1284*7? 


(272) 


(273) 


(274) 
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f  \ 

x„  -  X 


V  J 


RT 

F 


■c 


a 


k  , 


=  1.1284 4r 


(275) 


Note  that  the  left-side  of  Eq.  (275)  is  taken  to  be  the  measure  of  voltage  change  at  any  time 

t  =  .  It  is  suggested  that  the  result,  Eq.  (275)  should  be  used  only  for  t  values  such  that 

a 


8  4\[Dt  (7)7  <— 

Y  =  Ym  =  —  =  — - =  4.1—  =  4v  r  <  0. 1  to  satisfy  the  condition  of  short  penetration  of  Li+  (as 

a  a  V  a~ 

Li+  +  e  )  in  the  cathode  active  material  spherical  particle  from  the  electrolyte-active  material 
particle  interface.  We  may  express  Eq.  (275)  as 


f  \ 

x,  -  x 


V  J 


RT 
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K, 
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Q  v 

_ 1 

=  1.1284n/7 


(276) 


where  (Al/,))(  is  indicative  of  the  voltage  change  associated  with  the  lithium  ion  diffusion 
inside  the  cathode  active  material. 


From  Eq.  (256),  one  can  obtain 


x„  =  x  + 


yf; T 


yfr  ( Ja ) 
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V  P cathode^  J 


=  x. + 


2.1— 
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v  j 


M 


V  F cathode  ^  J 


(277) 


Equation  (277)  may  also  be  used  to  show  the  effect  of  transient  accumulation  of  lithium  inside 
the  cathode  active  material  particles  at  a  fixed  current  density  |E  |  on  the  reversible  voltage  of 
the  cell  through  the  application  of  Eq.  (266)  or  (143). 

4.1.9  Further  Treatment  of  the  Reversible  Potential  of  the  Lithium-Based  Cell,  Eq.  (144) 

From  Eq.  (144),  one  can  obtain  the  following  equation: 


E  =E 

rev  rev 


RT  (  a\  RT  f 

+^ln(x 


In 
F  l  L 


Yu-, 


in  MvPc 


+  lnx 


(278) 


After  entering  the  cathode  active  material  through  the  electrolyte-active  material  particle 
interface,  lithium  ion  moves,  with  its  “captured  electron”  present  in  its  near  vicinity.  This 
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situation  may  be  represented  as:  Li+  +  e  □  Li  inside  the  cathode  active  material  acting  as  the 
solid-state  solvent.  Then,  the  following  expression,  similar  to  that  for  a  liquid-phase  electrolytic 
solution  [10],  is  written  as  an  approximation  to  Inf yUAn M  Pc 


r  n 

z+z_asjl 

^Li-in  My  Pc 

\  +  Ba4l 

+  pxx 

II 


I 


(279) 


Term  I  of  the  right-hand  side  of  Eq.  (279)  is  assumed  to  account  for  the  lithium  ion  interactions 
with  the  surrounding  charges  through  the  coulomb  type  (long-range)  forces  and  the  second  term 
II  is  assumed  to  account  for  the  short-range  interactions  of  lithium  ion  with  the  surrounding 
charges. 


Term(I)  = 


z.z. 


a\l  (-1  )(-l  )«\//  asjl 


l  +  Basfl  \- Basil  \~  Basil 


(280) 


2  p  J^P  ^ 

where  a  =  lithium-ion  diameter;  parameter,  B  =  ,  — —  ■  F;  a  =  — — — - — 


and 


sRT  8  xL(sRT) 

I  =  —  ^  zfmi  =  ~[^Li+mLi+  +  zl  m,  ]  =  —  [(l)m  +  (l)m]  =  m  =  molality  of  the  intercalated  lithium 


in  the  solid  state  cathode  material  as  solvent  = 


M 


.  L  =Avogadro’s  number  -  6.0223x10 
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perg-mole,  F  =Faraday’s  constant  =  96487  coulomb  per  g-equivalent,  R  -  gas  constant 
=  8.3 14  J  /  mol  ■  K ,  T  —  temperature  [  AT] ,  t:  =permittivity  of  the  solid-state  cathode  active 

material  =  Dso  -  (dielectric  constant  of  the  solid-state  cathode  active  material)  (permittivity  of 


free  space),  so  =8.8542x10 


,4  Farad  0  0.l0  ...  I4  coulomb  volt  1 


=  8.8542x10“ 


cm 


cm 


s,  -mole  of  the  mobile  or  intercalated  lithium  ions 

x  =  ^ - ,  M 

g-mole  of  the  cathode  active  material 

material,  po  =  density  of  the  cathode  active  material. 


cathode 


-  mol.  wt.  of  the  cathode  active 


Using  the  information  in  Eq.  (280),  Eq.  (279)  is  expressed  as: 
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(281) 
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Inserting  for  In  7Li.inMPc  into  Eq.  (278)  leads  to: 


Equation  (282)  is  now  expressed  as: 


v  RT  , 

Erer+—  ln  *  ~  ‘ 

F 


RT  (  /?T  \ 

=  (/)  -  yA(  (283) 


This  result  suggests  that  one  should  experimentally  determine  Erev  vs.  x,  i.e., open-circuit 
voltage  vs.  x  data  at  the  cathode  active  material  particle  equilibrium  state  at  each  x  value  at  a 
fixed  temperature.  Then,  one  should  plot 


r  ^RT  ^ 

Erer+—  ln  *  “  ‘ 

F 


vs.  x  on  a  linear  graph  paper.  If  the  plot  is  predicted 


to  be  a  straight  line  with  slope  =  — — J/?,  ’  1'lcn’  A  's  determined  for  its  application  in  Eq. 

(281),  and,  hence,  in  Eq.  (278)  for  future  computations  of  Erev  vs.  x  at  the  temperature  for 
which  /?,  has  been  estimated.  Note  that  x°  is  the  x  value  to  be  selected  corresponding  to  a 
E°rev  value  and  also  x  should  be  greater  than  x°  value. 


4.1.10  Fractional  Utilization  of  a  Cathode  Active-Material  Particle  as  a  Function  of  Time 
at  a  Constant  Electrolyte- Active  Material  Interfacial  Current  Density 

The  fractional  utilization  of  the  cathode  active  material,  U,  is  defined  as  a  function  of  time  as 
follows. 
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^amount  of  lithium  transferred  into  a  cathode  active  material  particle  of 
radius,  a,  in  time,  t,  during  the  cell  discharge  period  at  a  constant  lithium- 
Ip  on  flux  or  current  density  at  the  electrolyte-active  material  interface 
^amount  of  intercalated  lithium  the  active  material  ' 
vparticle  of  radius,  a,  can  hold  up  to  the  saturation  level  J 
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where  Js  =  lithium  ion  flux  and  if  =  constant  electrolyte  active  material  interfacial  current 
density. 


Note  that  V  should  be  <  tssa,d ,  the  time  in  which  the  lithium  concentration  at  the  particle 
surface  reaches  the  saturation  level,  i.e.,  ( xs  -  xini  )  becomes  (xsald  -  .  This  t[aU>  can  be 

determined  using  Eq.  (169)  in  the  following  form. 
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(287) 


One  should  determine  rssa,d  by  solving  Eq.  (287)  for  a  fixed  value  of  if .  Then,  through  the 


definition  of  rf""/  =  — 1 


sat'd 


Q^''£Sat'd 

;  one  determine  Patd  =  — = — 
s  D 
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5.  CONCLUSIONS 


To  help  explain  AFRL/RZPS  in-house  experimental  results  on  lithium-based  cells  using 
dilithium  phthalocyanine  as  the  active  cathode  material,  Dr.  Sarwan  Sandhu  from  the  University 
of  Dayton  was  tasked  to  develop  a  model  formulation  and  simulation  for  these  types  of 
electrochemical  cells.  The  experimental  and  modeled  cells  consisted  of  lithium  foil  anodes, 
lithium  tetrafluoroborate  in  poly(ethylene  oxide)  electrolytes,  and  metal  phthalocyanine  based 
cathodes.  The  developed  formulation  covered  the  various  subtopics  relating  to  the  performance 
and  processes  occurring  at  the  cell  cathode,  electrolyte  separator,  and  cell  anode. 

Determination  of  the  transient  cathode  performance  under  discharge  was  treated  extensively  due 
to  the  assumed  solid-state  diffusion-limited  cell  behavior.  Boundary  conditions  of  a  constant 
concentration  or  current  density  on  the  cathode  active  material  particles  were  used.  Solutions 
consisted  of  the  cathode  performance  for  a  given  surface  concentration  or  for  a  given  cell  current 
density. 

The  reversible  cathode  potential  was  determined  as  a  function  of  the  amount  of  lithium 
intercalated  into  the  phthalocyanine  material. 

Equations  for  the  electrode  kinetics  at  the  cell  cathode  and  anode  were  developed  to  help 
determine  the  cathode  overpotentials  or  activation  polarization  losses  occurring  during  cell 
discharge. 

Further  effort  is  needed  on  using  the  developed  formulation  and  applying  it  to  existing  and  future 
experimental  results  on  lithium-based  phthalocyanine  solid-state  cells. 
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